Magnetic States at Short Distances 

Horace W. Crater^* | and Cheuk-Yin Wong^'^'f 

'^ The University of Tennessee Space Institute, Tullahoma, Tennessee 37388 

^Department of Physics and Astronomy, University of Tennessee, Knoxville, TN 37996 

'^Physics Division, Oak Ridge National Laboratory, Oak Ridge, TN 37831 

(Dated: February 24, 2013) 

The magnetic interactions between a fermion and an antifermion of opposite electric or color 
charges in the ^Sq^ and '^P^^ states with J = are very attractive and singular near the origin 
— . 1 ■ and may allow the formation of new bound and resonance states at short distances. In the two 

I , body Dirac equations formulated in constraint dynamics, the short-distance attraction for these 

^— ^ ■ states for point particles leads to a quasipotential that behaves near the origin as —a^/r'^, where 

^vj ' a is the coupling constant. Representing this quasipotential at short distances as A(A + l)/r^ with 

' A = (—1 + y/l — 4a^)/2, both ^S^^ and '^Pq'^ states admit two types of eigenstates with drastically 

f^ . different behaviors for the radial wave function u = rip. One type of states, with it growing as 

^M • ^-^+1 at small r, will be called usual states. The other type of states with u growing as r^"^ will be 

^^ ' called peculiar states. Both of the usual and peculiar eigenstates have admissible behaviors at short 

f~^ , distances. Remarkably, the solutions for both sets of ^ So states can be written out analytically. The 

fT^ ' usual bound ^So states possess attributes the same as those one usually encounters in QED and 

QCD, with bound state energies explicitly agreeing with the standard perturbative results through 
order a'*. In contrast, the peculiar bound ^So states, yet to be observed, not only have different 
behaviors at the origin, but also distinctly different bound state properties (and scattering phase 
*Y^' shifts). For the peculiar ^ So ground state of fermion-antifermion pair with fermion rest mass m, 

O „ the root-mean-square radius is approximately 1/m, binding energy is approximately (2 — y/2)m, 

^ . and rest mass approximately \f2m. On the other hand, the (n-f 1)^50 peculiar state with principal 

<~i ' quantum number (n -I- 1) is nearly degenerate in energy and approximately equal in size with the 

n}So usual states. For the ^Po states, the usual solutions lead to the standard bound state energies 
and no resonance, but resonances have been found for the peculiar states whose energies depend on 
the description of the internal structure of the charges, the mass of the constituent, and the coupling 
constant. The existence of both usual and peculiar eigenstates in the same system leads to the non- 
self-adjoint property of the mass operator and two non-orthogonal complete sets. As both sets of 
\f^ ' states are physically admissible, the mass operator can be made self-adjoint with a single complete 

^— s ' set of admissible states by introducing a new peculiarity quantuin number and an enlarged Hilbert 

• , space that contains both the usual and peculiar states in different peculiarity sectors. Whether 

C^ ■ or not these newly-uncovered quantum-mechanically acceptable peculiar ^So bound states and '^Po 

>— ' ' resonances for point fermion-antifermion systems correspond to physical states remains to be further 

investigated. 
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I. INTRODUCTION 

It is well known that for some combinations of the spin configurations and orbital motion the magnetic interaction 
can be strongly attractive and singulaiQ at short distances [2|-13| • We can illustrate this by a classical example as shown 
schematically in Fig. 1(a) where a positive charge g"'" is making a circular orbit about a fixed negative charge q^ 
whose spin s{q~) is pointing in a direction opposite to the orbital angular momentum of g"*" [J|. In the external field 
problem, (e.g., Fermi's treatment of hyperfine structure), the charged particle q~ at rest with a magnetic moment 
n{q~) generates a vector potential A = fJ^iq^) x r/r^ which acts on the other particle, q'^ . Such a "magnetic" 
interaction can be very attractive when the spins and the orbital angular momentum are oppositely aligned, as shown 
in the configuration of {q~^q~) in Fig. 1, where the vector potential A, arising from the q~ magnetic dipole moment 
fi{q^), is parallel to the g+ orbital momentum p. The interaction (— p • A) from q^ acting on q^ is attractive and is 
proportional to [L{q'^) ■s{q~)]/r^ that is quite singular in nature. At short distances it may overwhelm the centrifugal 
barrier that is proportional to l/r^. Similarly, the interaction from q^ acting on q~ will be likewise attractive and 
singular if the spin of the s(q+) is parallel to the electron spin s{q—) and pointing in the same direction, resulting 
in the total spin of the q'^q^ system aligning opposite to the orbital angular momentum, as in the '^Pq state with 
S- = 1, L = 1, J = 0, P = +1, and C = +1. 

The '^Pq'^ state is not the only state with a strong magnetic interaction. One can envisage classically another spin 
configuration, the ^Sq~^ state, that also has attractive and singular magnetic interactions. As illustrated schematically 
in Fig. 1(b), a fermion q~ with an electric or color charge interacts with an antifermion q~^ of opposite electric or 
color charge with spins s{q~) and s((7"'") pointing in a opposite directions in the ^Sq state configuration. With the 
spins opposite to each other, the magnetic moments of q~ and g"*" are parallel to each other. The interaction between 
the magnetic moments is @ i^int = ~(87i'/3)/.tq- • fi^j+S^r), which is attractive and singular at short distances. The 
strong and singular magnetic interaction may overcome other repulsive interactions and may allow the formation of 
bound states of the fermion and antifermion system at short distances. For brevity of notation, the quantum numbers 
P and C in and ^Sq^ and ^P^^^ will be understood. 
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FIG. 1. (a) The schematic picture of the ^Pq state spin configuration and the orbital motion of a negative charge q~ and a 
positive charge g"*" that can lead to a strong magnetic attraction at short distances. Here, /x(g ) is tire magnetic moment of 
the charge q^ arising from its spin s{q^). (b) Tire scliematic picture of the spin configurations of q^ and g+ in the ^^o state. 

Previously, one of us (CYW), in collaboration with R. L. Becker, studied the (e+e~) system using the Kemmcr- 
Fermi-Yang equation [y] with interactions consisting of the Coulomb interaction and the vector (magnetic) interaction. 
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■^j in connection with a possible scalar '^Pq magnetic resonance [J]. The interest was to 
investigate whether there could be a resonance at the mass of 1.579 MeV that might explain the anomalous positron 
peak in heavy- ion collisions near the Coulomb barrier [7[ . The experimental evidence for the anomalous positron peak 
later turned out to be negative when greater statistics were accumulated [g|. Nevertheless, it remains of interest to 
study the behavior of the two-body system at short distances and see how the attractive magnetic interaction in the 
^Pq state may reveal itself in some observable properties. 

While the use of the Kcmmer-Fcrmi-Yang equation with a two-body magnetic interaction is useful to motivate an 
approximate description [J] , a consistent relativistic description of the two-body interaction at short distances can be 
found in the relativistic two body Dirac equations (TBDE) formulated in Dirac's constraint dynamics [9l-[l3|. These 



A potential is quantum mechanically singular if it is more attractive than 
See Q. 



-l/4r^ at the origin in the context of {~d?/dr'^ — l/4r-^). 



rclativistic two body Dirac equations give a good description to the entire meson mass spectrum (excluding most 
flavor-mixed mesons) with constituent world-scalar and vector potentials depending on just two or three invariant 
functions, in previous relativistic quark- model calculations [13l4l5| . 

The application of the TBDE equations to two-body bound and resonance states in quantum electrodynamics 
has intrinsic merits. In Ref. [1^1, the properties of these TBDE equations that made them work so well for the 
rclativistic quark model were investigated by solving them nonperturbatively (i.e. analytically or numerically) in 
quantum electrodynamics (QED), where order a^ perturbative solutions are well known. The two coupled Dirac 
equations in the constraint formalism depend on Lorentz-covariant potentials between the two constituents and act on 
a 16-component wave function. An exact Pauli reduction led to a second-order relativistic Schrodinger-like equation 
for a reduced four-component wave function with an effective interaction containing all the dependencies on spin, 
orbital angular momentum, and tensor operators. We were able to solve the TBDE nonperturbatively (analytically 
or numerically) as well as pert urbatively because the spin dependent short-distance components of the effective 
interaction are not singular J12l Il6j . The situation is very different from the approximate Fermi-Breit forms, which 
contain singular potentials and necessitate the introduction of arbitrary short-distance cut-off parameters. The spin 
dependence of the relativistic potentials in the exact Schrodinger-like equation arises naturally from the relativistic 
reduction procedure and it incorporates detailed minimal interaction and dynamical recoil effects, characteristic of 
field theory. We shall also use the term "quasipotential" to represent this effective, non-singular interaction. 

To obtain the interaction used in the TBDE formalism, we first determined the relativistic quasipotential to the 
lowest order in a for the Schr odinger-like equation in Ref. [l6| by comparing the effective interaction with the 
interaction derived from the Bethe-Salpeter equation. This, in turn led to an invariant Coulomb-like potential A{r) = 
~a/r, where a is the coupling constant. Insertion of this information into the minimal interaction structures of the 
two body Dirac equations then completely determined all aspects (spin-dependent as well as spin-independent) of the 
interaction. (In |17| we gave a procedure to construct the full 16-component solution to our coupled first-order Dirac 
equations from a solution of the second-order equation for the reduced wave function.) 

Next, we showed that both the quantum mechanical perturbative and the TBDE non-perturbative treatments 
(i.e. analytic or numerical) yield the standard spectral results for QED and related interactions through order a^ 
. Such an agreement depends crucially on the inclusion of the coupling between various components of our 16- 
component Dirac wave functions and on the short-distance behavior of the relativistic quasipotential in the associated 
Schrodinger-like equation. We then examined the speculations [J| whether the quasipotentials (including the angular 
momentum barrier) for some states in the e+e" system may become attractive enough at short distances to yield a 
pure QED resonance corresponding to the anomalous positron peaks in heavy- ion collisions 0. For the ^Pq state 
we found that, even though the quasipotential becomes attractive and overwhelms the centrifugal barrier at short 
distances, the spatial extension of the attraction is not large enough to hold a resonance at the energy of 1.579 MeV 
\Wlj . This result contradicted predictions of such states by other authors [T^l based on numerical solutions of three- 
dimensional truncations of the Bethc-Salpeter equation, for which the entire QED bound state spectrum has been 
treated successfully through order a^ only by perturbation theory. 

In this paper we return to this problem of the magnetic resonance and magnetic states, not motivated so much 
by new experimental data as by a discovery of an additional peculiar solution of the TBDE overlooked in the earlier 
work in Ref. [l6|. Our examination of the two body Dirac equations reveals that at short distance for both ^5*0 and 
'^Pq states, the magnetic interactions is indeed quite strong. As a consequence, they counterbalance other repulsive 
interactions to result in a quasipotential for these states that behaves as —a^jr^ at short distances. 

In standard quantum mechanics for central interactions including the angular momentum barrier L(L + l)/r^ for 
states with L ^ at short distances, one generally retains only one of the two solutions for the radial part of the wave 
function, u = ri/j, the one that grows with distance as (~ r'"^^), dropping the other solution(~^ r^^) as being too 
singular. If we likewise represent the quasipotential as A(A + 1)/?'^ with A = (—1 + ■\/(l — 4a^)/2, it leads to a short- 
distance solution that behaves as r^+^, which we call the usual solution, in addition to a solution, whose radial part 
grows as r"'*', which we call the peculiar solution. However, both usual and peculiar states have quantum-mechanically 
acceptable behaviors at short distances, as the wave functions at short-distances are square-integrable. 

In the case of the spin singlet ^^o states, the eigenstates and eigenenergies can be obtained analytically and are 
found to encompass both usual and peculiar states. We find usual bound states with attributes the same as those one 
usually encounters in QED and QCD, explicitly agreeing with the standard perturbative results through order a*. In 
contrast, the peculiar ^5*0 ground state of a fermion-antifermion pair with a fermion rest mass m has a root-mean- 
square radius approximately l/m, a binding energies approximately Bp^{2 — \/2)m, and a rest mass approximately 
\/2m. However, the (n + 1) th ^5*0 peculiar state is nearly degenerate in energy and approximately equal in size with 
the nth usual ^Sq. 

The existence of both usual and peculiar eigenstates in the same system brings with them conceptual and mathemat- 
ical problems of the non-self-adjoint property of the mass operator and the over-completeness of the set of eigenstates. 
We resolve these problems by the introduction of a new quantum number, the peculiarity quantum number, that 



makes the mass operator self-adjoint and the combined set of usual and peculiar states a complete set in an enlarged 
Hilbert space. 

In the case of the ^Pq states, both of the usual and peculiar solutions reflect the overwhelmed centrifugal barrier 
and so differ substantially from the r^+^ and r~^ behaviors at short distances respectively. As a peculiar state radial 

,2 

wave function u rises from the zero value at the origin as r~ ~ r" , the strongly attractive magnetic interaction has 
the tendency of bending the wave function in such a way to allow for the possibility of a resonance. Furthermore, 
as the quasipotential obtained through the relativistic reduction is sensitively energy dependent, we can explore the 
behavior of the two-body system over a larger domain of energies. We find that the usual solutions lead to no resonant 
behavior, but the peculiar solution can lead to a ^Pq resonance whose phase shift changes by tt at an appropriate 
energy, depending on the description of the internal structure of the charges, the mass of the constituent, and the 
coupling constant. 

This paper is organized as follows. In Sec. II we give a review of the two body Dirac equations of constraint 
dynamics. For those readers who are already familiar with the constraint approach we refer them to the TBDE 
given in Eq. P^ and their Schrodinger-like Pauli reduction given in Ec^. ([T7)) . We specialize to electromagnetic-like 
interactions only in this paper. We give in Sec. Ill the single-component radial forms of Eq. (|17p relevant to this paper. 
In Sec. IV we examine both solutions for the ^5*0 states. In addition to examining new bound state solutions, we show 
how the -^Sq wave functions for positive energies (and their corresponding phase shifts) can be determined analytically 
in terms of Coulomb wave functions for noninteger angular momentum. This is done for both the usual and peculiar 
solutions. We explain why and how we introduce of a new quantum number, which we call the peculiarity quantum 
number, to solve the problems of the non-self-adjoint property of the mass operator and the over-completeness of 
the set of eigenstates. In Sections V we examine the short distance behaviors for the ^Pn state for the usual and 
peculiar solutions. In Sec. VI we discuss the variable phase shift formalism of Calogero [l^ and outline how we use 
it for our phase shift analysis. Since the short distance behavior of the '^Pq quasipotential is the same as that of 
^5*0 quasipotential, we can use those same ^5*0 Coulomb wave functions as reference wave functions in that region 
to compute phase shifts. There is, however, an additional term (proportional to S{r)) that does not appear in the 
extreme short distance region for the ^Sq quasipotential. Even though this term does not contribute in the case of 
the phase shift for the usual solution, its contribution to the phase shift calculations for the peculiar solution must be 
considered. In Sec. VI we discuss our numerical results and in Sec. VII our conclusions. Various technical results are 
presented in the appendices. In Appendix A we give an outline of the details on the relation between the two-body 
Dirac equations and their Pauli reduced Schrodinger forms. In Appendix B we present the radial forms of those 
Schrodinger-like equations for a general angular momentum coupling. In Appendix C we present details of the ^Sq 
usual and peculiar bound states. In Appendix D we review the connections between the Coulomb wave functions for 
noninteger angular momentum index. Appendix E presents a review of the variable phase method of Calogero [l9[ 
for our problem. 



II. TWO BODY DIRAC EQUATIONS 

We briefly review the two body Dirac equations of constraint dynamics [lll4l5l |20| providing a covariant three 
dimensional truncation of the Bethe Salpeter equation for the two body system. Sazdjian [23 - l24l | has shown that the 
Bethe-Salpeter equation can be algebraically transformed into two independent equations. The first yields a covariant 
three dimensional eigenvalue equation which for spinless particles takes the form 

Hlo+H20+2$W(xi,.X2)=0, (1) 

where Hio = pf + ml . The quasipotential $ is a modified geometric series in the Bethe-Salpeter kernel K such that 
in lowest order in K 

$ = Tii5{P ■ p)K, (2) 

where P = pi + P2 is the total momentum, p = /i2Pi — /^iP2 is the relative momentum, w is the invariant total center 
of momentum (cm.) energy with P" = —iv. The fii must be chosen so that the relative coordinate a; = xi — a;2 and 
p are canonically conjugate, i.e. /ii + /i2 = 1- The second equation, Eq. ([2]), overcomes the difficulty of treating the 
relative time in the center of momentum system by setting an invariant condition on the relative momentum p, 

{nio-n2o)'^{xi,X2) = = 2P •p*(xi,X2). (3) 



Note that this impUes p'"* = p^* = (77'''' + P''P'^)pi,'9 in which P^ = P^/u; is a time hkc unit vector (P^ = -1) in 
the direction of the total momenturro. 

One can further combine the sum and the difference of Eqs. ([l} and ((3)) to obtain a set of two rclativistic equations 
one for each particle with each equation specifying two generalized mass-shell constraints 

H.*(xi,.T2) = (p?+m2 + $)^(a:;i,a;2) = 0, ^ = l,2, (4) 

including the interaction with the other particle. These constraint equations are just those of Dirac's Hamiltonian 
constraint dynamics for spinless particles [^, [l3, l2l|. In order for Eq. ([4]) to have consistent solutions, Dirac's 
constraint dynamics stipulate that these two constraints must be compatible among themselves, [Hi,H2]'i ~ 0, that 
is, they must be first class constraints. This requires that the quasipotential $ satisfy [pf —P2, ^]'^ = 0. Working out 
the commutator shows that for this to be true in general, $ must depend on the relative coordinate a: = xi — a;2 only 
through its component, x±, perpendicular to P, 

x>i^{7J^''' + P^Pn{xl~X2).. (5) 

The invariant x^ = r^ becomes r^ in the cm. frame. Since the total momentum is conserved, the single com pon ent 
wave function ^ in coordinate space is a product of a plane wave eigenstate of P and an internal part '4'{^±) 25^1. 



We find a plausible structure for the quasipotential $ by observing that the one-body Klein-Gordon equation 
(p^ -I- TO^)-0 = (p^ — e'^ + m?)'ip ~ takes the form (p^ ~ e'^ + m? + 21118 + S"^ + 2eA — A'^)il' ~ when one introduces 
a scalar interaction and time-like vector interaction via the minimum substitutions m ^ m -\- S and e — > £ — A. 
In the two-body case, separate classical [26[ and quantum field theory |24| arguments show that when one includes 
world scalar and vector interactions then $ depends on two underlying invariant functions S{r) and A{r) {r = \/x\) 
through the two body Klein-Gordon-like potential form with the same general structure, that is 

^ = 2m.^S + S'^ + 2e.ujA-A^. (6) 

Those field theories further yield the cm. energy dependent forms 

m,u, = Tnim2/w, (7) 

and 

£w = (w^ -m\- m\)/2w, (8) 

ones that Tododov [l^, [231 introduced as the rclativistic reduced mass and effective particle energy for the two-body 
system. Similar to what happens in the nonrelativistic two-body problem, in the rclativistic case we have the motion of 
this effective particle taking place as if it were in an external field (here generated by S and ^4) . The two kinematical 
variables ([7]) and ([8]) are related to one another by the Einstein condition 

el-ml = h\w), (9) 

where the invariant 

b^{w) = [w^ - 2w^{m\ + ml) + {m\ - mlf)/Av?, (10) 

is the cm. value of the square of the relative momentum expressed as a function of w. One also has 

h\w)^el-ml=el-ml (11) 

in which £1 and £2 are the invariant cm. energies of the individual particles satisfying 

£1 -I- £2 = w, £1 - £2 = [rni - mD/w. (12) 

In terms of these invariants, the relative momentum appearing in Eq. ^ and ^ is given by 

p^ = {e2Pt - e,p^^)lw, (13) 



^ We use the metric -q^^ = •tf''^ = rf''^ = — r;"" = 1. 

^ We use the same symbol P for the eigenvalue so that the w dependence of m^ and e-w in Eq. ^ is regarded as an eigenvalue 
dependence. The wave function \l/ can be viewed either as a rclativistic 2-body wave function (similar in interpretation to the Dirac 
wave function) or, if a close connection to field theory is required, related directly to the Bethe Salpoter wave function x t>y |23ll 
* = ~TTiS{P ■ p) Hiox = -TTiSiP ■ p)'H20X- 



so that 1^1+ fJ,2 = (si + £2)/^ = 1. In [2^ the forms for these two-body effective kinematic variables are give n sound 
justifications based solely on relativistic kinematics, supplementing the dynamical arguments of [26| and |24l |. 

This covariant and useful three-dimensional truncation of the Bethe-Salpeter equation has been extended to the 
case of a two-fermion system where the two constraint equations become the two body Dirac equations (TBDE) 

[iMl 

Sii' = 751(71 • {Pi - Ai) + mi + 5i)* = 0, 

52V' = 752(72 ■ {P2 - A2) +1712+ ^2)* = 0. (14) 

Here ^ is a sixteen component wave function consisting of an external plane wave part that is an eigenstate of P and 
an internal part ip = 4'{x±)- The vector potential A'' was taken to be an electromagnetic-like four- vector potential 
with the time-like and space- like portions both arising from a single invariant function A{r). |f| The tilde on these 
four-vector potentials indicate that they are not only position dependent but also spin dependent by way of the 
gamma matrices. The operators Si and ^2 must commute or at the very least [Si,S2]'ip = since they operate on 
the same wave functioiij. This compatibility condition gives restrictions on the spin dependencies of the vector and 
scalar potentials, 

if = if(A(r),pi,P,«^,7i,72), (15) 



in addition to requiring that they depend on the invariant separation r = \l x\ through the invariant A(r). The 
covariant constraint (jS]) can also be shown to follow from Eq. ((T4)) . We give the explicit connections between A^ and 
the invariant A(r) in Appendix A. (A similar dependence occurs for Si on S{r).) The general structural dependence 
on A{r) and S{r) and the spin dependence of Af , St is a consequence of the compatibility condition [Si,S2]4' = 0. 

The Pauli reduction of these coupled Dirac equations lead to a covariant Schr odinger-like equation for the relative 
motion with an explicit spin-dependent potential $, 



Pi + $(5(r), A(r),px, P, w, ai, (72) 1 V+ = b'iw)ij+, (16) 

with b^ (w) playing the role of the eigenvaluqj. This eigenvalue equation can then be solved for the four-component 
effective particle spinor wave function "0+ related to the sixteen component spinor ip(x±) in appendix A. 

In Appendix A we outline the steps needed to obtain the explicit cm. form of Eq. ([T6| . That form is [23], j30l |. 

[iM3 

{ p^ + $(r,mi, 7712, w,cri, 0-2) }V'+ 

= {p2 + 2m.^5 + S^ + 2eu,A -A^ + <^d 

+ L ■ (0-^+cr2)$S0 + (Ti-f(72-fL ■ {(7^+(72)^SOT 
-|-(Ti-cr2*&SS + (3ci-fcr2 ■ f — (Ti-(J2)'^T 
+L ■ {(^i-(^2)^sOD + *-^ • CTi Xcr2'5sox}'0+ 

= 6^^+, (17) 

where the detailed forms of the separate quasipotentials $i are given in Appendix A. The subscripts of most of the 
quasipotentials are self explanatorjQ. After the eigenvalue iP' of ([T7| is obtained, the invariant mass of the composite 
two-body system w can then be obtained by inverting Eq. (jlOp . It is given explicitly by 



b'^+ml + ^b^+rnl. (18) 

For this reason we call the operator that appears to the left of Eq. (flT)) the invariant mass operator. The structure 
of the linear and quadratic terms in Eq. ()17p as well as the Darwin and spin-orbit terms, are plausible in light of 
the discussion given above Eq. (j6]), and in light of the static limit Dirac structures that come about from the Pauli 
reduction of the Dirac equation. Their appearance as well as that of the remaining spin structures are direct outcomes 
of the Pauli reductions of the simultaneous TBDE Eq. ([T4| . In this paper we take the scalar interaction S{r) = 0. 



In particular, in a porturbative context that would mean that these aspects of A'^ were regarded as arising from a Feynman gauge 
vertex coupling of a form proportional to 7^72^ A . 

The 75 matrices for each of the two particles are designated by 75; i = 1,2. The reason for putting these matrices out front of the 
whole expression is that including them facilitates the proof of the compatibility condition, see Illl. l25ll . 

Due to the dependence of <& on u), this is a nonlinear eigenvalue equation. 

The subscript on quasipotential ^0 refers to Darwin. It consist of what are called Darwin terms, those that are the two-body analogue 
of terms that accompany the spin-orbit term in the one-body Pauli reduction of the ordinary one-body Dirac equation, and ones 
related by canonical transformations to Darwin interactions |26| . 13711 , momentum dependent terms arising from retardation effects. The 
subscripts on the other quasipotentials refer respectively to SO (spin-orbit), SOD (spin-orbit difference), SOX (spin-orbit cross terms), 
SS (spin-spin), T (tensor), SOT (spin-orbit-tensor). 



III. TBDE SINGLE COMPONENT WAVE EQUATIONS. 

The 4 component two-body wave function 1/'+ of the above PauU-form ( [17)) of the TBDE can be conveniently 
represented by spin-singlet 5 = and spin-triplet 5 = 1 components with quantum numbers {J, L, S} and basis wave 
functions 

{r\wJLS) ^ ^jLs{r) = ^^^^^^YjM{i). (19) 

r 

In general, the singlet and triplet states arc coupled. However, we see from Appendix B that for the case of equal 
masses and certain angular momentum states, the spin singlet and spin triplet components decouple, and the TBDE 
reduce to a single component equation. 

Specifically, for the spin-singlet 5 = state with J = L, (the ^Jj state), the TBDE is 

{~;^ + "^ ^^1 + 2£"'^ - ^' + *i5-3$ss}m,/,/o = b^ujjo, (20) 

where, using the results in Appendix B, the magnetic interaction — 3$55 is 

1 3 / A' \^,, 1 



- 3^ss = -3$ss(A A', v^A) = -3(- - - [^;;^))ii^r=^jj^ + V^^^AJ^) - 2) 

'2? \w-2a) ^yT~~2l7^ /^i-y \w-2a) ^ W-2A ^ ^' 

= --^d{A,A',V^A), 

which is attractive and singular, as we discussed in the Introduction. At large distances and for A — —a/r potential, 
'V^A = 4:'KaS{r) and the spin-spin interaction indeed becomes a singular interaction as described in Q. In addition 
to the magnetic spin-spin interaction, there is also the repulsive Darwin quasipotential $£>. In the ^Jj state, the 
attractive magnetic spin-spin quasipotential in the spin-singlet configuration exactly cancels the repulsive Darwin 
quasipotential, 

-3$ss + $D = 0. (22) 

As a result of this remarkable cancellation, the eigenvalue equation for the ^ Jj state in Eq. (PT|) becomes simply 

{-:r2 + "^2 + 2£«-^ - A^}uj.ja = b^ujjo. (23) 

ar'' r^ 

Of all spin-singlet states, only in the ^Sq states (J = L = 0) do the effects of the quasipotential and the absence of 
a centrifugal barrier make the combined quasipotential strongly attractive at short-distances. This, of course, would 
not happen were it not for the highly attractive spin-spin interaction discussed in the Introduction and in Eq. (|2ip . 
Among the spin-singlet states with different J quantum numbers, we shall therefore focus our attention only on the 
^5o states. 

For the spin-triplet 5=1 states, there are two states with single component radial equations. The first is the ^Jj 
state whose radial equation takes the form (J > 1) 

, d^ J(J+1) „ , ,2 2A' 0/ A' Y V^A , ,2 



dr^ r2 •" r{w-2A) \w~2Aj w - 2A 

The second is the '^Pq equation which takes the form 

r d"^ 2 , ,, 8A' f A' Y 2V^A , , , , 

{-T^ + — + 2£»^-^ -^ ^7T+^ ^ + ^}"oii = &'woii. 25 

dr'^ r^ r [w — 2A) \w — 2AJ w — 2A 

Of the two spin-triplet cases, only in the ^Pq states ( J = 0, L = 1) do the combined effects of the quasipotentials 
become so strongly attractive at short-distances that they overwhelm the presence of the centrifugal barrier. As 
discussed in the Introduction, this is due to the highly attractive spin-orbit interaction ( "magnetic" interaction) when 
the total spin and the orbital angular momentum are oppositely aligned. In that case, the competing effects of both 
the short-distance attraction and the presence of the potential barrier raise the question whether the attraction is 



strong enough to hold a resonance state in the continuum. Among the spin-triplet states with different J and L 
quantum numbers, we shall therefore focus our attention only on the '^Pq states. 

In the last term of the quasipotential in Eq. (|25p . the quantity V'^A is related to the particle charge density, p(r), 
seen by each of the two particles by 

VM(r) = 47rap(r). (26) 

Therefore, the equation for the two-body relative wave function for the ^Pq state becomes 

(P 2 ^ ,„ 8A' f A' \^ 87rap(r)1 ,, , , 

-^^^v^^ 2^-^ - ^ - 7F^2i) + « [-^;r^A) + -^ir^V^^^ = ^ "-- ^^'^ 

As we shall see in this case, the attractive magnetic interaction overwhelms the centrifugal barrier, allowing the wave 
function to reach the short-distance region where the particle charge density p{r)^ if any, can be exposed for scrutiny. 
This is in contrast to the situation for states in which the centrifugal barrier dominates the short-distance region. 
In that case, the centrifugal barrier will prevent the wave function from reaching the short-distance region and the 
particle charge density will not make as a significant difference in observable quantified. We obtain the important 
result that the ^Pq quasipotential depends explicitly on the particle charge density p{r) at short distances. As a 
consequence, some observable quantities may depend more critically on the nature of the particle charge distribution 
and the forces binding the charge elements together. 

For the ^Pq state, it is convenient to separate out the centrifugal barrier 2/r^ and the quasipotential $ to write the 
above equation as 

<f 2 1 

^ + — + *(^)r"oii-=^'^woii, (28) 

where 



{- 



^. X „ . .2 8A' / A' Y Snapir) , ^ 

<^{r) = 2e^A- A^ ^rv + 8 ^—A + -irr- 29) 

^ ^ r[w-2A) \w-2AJ w - 2A ^ ' 

In our early work [Tg, we limited our attention to energy regions around 1.579 MeV for the (e+e") system and 
searched for ^Pq resonances whose wave functions start from the origin in the usual way. We found no resonance states. 
We return to this problem again including now an additional (peculiar) solution of the TBDE that was overlooked in 
the earlier work but has quantum-mechanically acceptable behaviors at short distances. 

In Eqs. ([271) ). both the gauge field A{r) and the gauge field source p{r) appear in the equation of motion for the 
wave function in the ^Pq state. The appearance of the fermion charge source distribution p{r) brings into focus the 
question whether it is sufficient to describe the magnetic interaction in the "^Pq state completely within quantum 
electrodynamics or quantum chromodynamics. Electrons in QED and quarks in QCD are taken to be point particles 
with no structure. It may be necessary to go beyond these field theories, to include additional auxiliary interactions 
that hold the charge elements together, in order to properly describe the internal structure of these particles. If these 
auxiliary interactions act on the charged elements of the fermion to hold them together, they can also act on the 
charged elements of the antifcrmion charge and will affect the ^Pq wave function in the interior region of the charge 
distribution p(r). 

The nature of these auxiliary forces holding the charged elements together is completely unknown, although there 
have been many attempts to carry out such an investigation. For example, in the Dirac's model of an electron, a 
surface tension from an unknown axillary interaction is invoked to hold the electric charged elements of an electron 
together [3ll - l3q . However, our knowledge on the internal structures of electrons and quarks remain very uncertain. 
We shall return to examine how such a lack of knowledge of the internal structures of these elementary quanta leads 
to uncertainties in the '^Pq magnetic resonance states in Sec. VC. 

IV. SOLUTIONS OF THE TWO BODY DIRAC EQUATIONS FOR THE ^5o STATE 

A. The ^So quasipotential 

We first consider the case of the ^Sq state of a point fermion-antifermion pair with electric or color charges interact 
through an electromagnetic-type interaction arising from the exchange of a single photon or giuon. The single photon 



Such would also be the case for ^^o states in which, due to the cancellation in Eq. I I22I I. the dependence on p{r) is only indirect or 
implicit through the altered form for A{r). 



annihilation diagram docs not contribute because the ^So state is a charge parity even state. We thus have 



a 



A = . (30) 

r 

For brevity of notation in this subsection, we shah abbreviate the radial wave function wojj=wooo a-s u. Equation 
(l20l) for u becomes 



d 2£tua a 



u = b-'u. (31) 



with a short distance (r << a/2ew) behavior given by 



-^^2-77^ = 0, (32) 



with solutions 



U+ r 





A = (-1 + v/r^^4^)/2, (33) 



^^^^ii±VT^4^)/2^ (34) 

both of which approach zero as r approaches zero. With these behaviors, the probability 



u2 



^l(fr = -^r^drdn = u\drdVl = r^^^^^^^'^drdfl, (35) 



r 



is finite for both signs. We call u^ the usual solution, and it behaves as 7-^+1 ~ r^ " for small a. We call u_ the 
peculiar solution, and it behaves as r~^ ~ r" for small a. Both of these behaviors are physically acceptable near the 
origin in the sense of (i) u(0)— > 0, and (ii) being square-integrable in the neighborhood of r ~ 0. We note that if 
the sign in front of a^ were positive or if we had non-zero angular momentum such that L{L + 1) — a^ > then the 
second or peculiar set of solutions are not physically admissible states. 

In |38| one finds a thorough discussion on the proper boundary condition for the radial wave function of the 
Schrodinger equation at the origin. They discuss several conditions that appear in the literature: (1) Continuity of 
R = u/r at r = 0, requiring u(0) = 0. (2) A finite differential probability in the spherical slice {r,r + dr), that 
is B?r^dr < cxd requiring u{t) -^ r'^'^^,s > —1 and again u(0) ~ 0. (3) Requiring a finite total probability inside 
a sphere of small radius a which allows a more singular behavior, namely u{r) — >■ f-i/^+e where e > is a small 
positive constant, which would also include a finite behavior of the norm. (4) Requiring time independence of the 
norm leading to u{r) — >■ cr~''+^, s < 1 which again leads to m(0) = 0. Reference [331 furthermore shows that the 
radial Schrodinger equation [—d'^/dr^ + l{l + 1)/^^ + 2mV{ry\u{r) = 2mEu{r) is compatible with the full Schrodinger 
equation (-V^ + 2mV{r))'^Yi^ = 2mE'^Yi„i if and only if the condition u(0) = is satisfied. This u(0) = 
condition is clearly satisfied for both solutions in Eq. ([M]). 

In SchifF's Quantum Mechanics [33, a solution similar to the peculiar one discussed here is examined for the case 
of the Klein-Gordon equation for the Coulomb system. He argues that what we call the peculiar solution can be 
discarded since the source of the Coulomb attraction is a finite sized nucleus of radius tq. In particular, he states 
that for r < vq for which the potential is finite all the way to the origin, matching at rg would rule out the peculiar 
solution. In our case, with point particles, the potential does not satisfy this condition. 

1. ^So Bound States 

The solutions of the ^Sq bound states can be obtained analytically. In Appendix C we show how we can obtain 
the two sets of ^Sq bound state solutions that correspond to the usual and peculiar short distance behaviors. The 
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respective sets of eigenvalues and normalized eigenfunctions for the state with total invariant cm. energy (mass) 'w±n 
and the principle quantum number n is 



w±n - my 2 + 2/^Jl + a^/{n ± ^1/4 - (^ - 1/2) 
u±„{r) = 



4et„,ar\ nJ. 



rij. J 2n[j.[n^ + X±y. „' \ n^' i 

where n'_j_ — rir + X± + 1 — n + X± and 

e„,^ = (w| - 2m'^)/2w±. (37) 

For the usual states u-|_„, the bound state eigenvalues w^n agree with standard QED perturbative results through 
order a^, 

w+n = 2m - ■ma'^/4n^ - ma'^/2n^{l - 11/32/n) + O(a^), n = l,2, 3, ... (38) 

For the set of peculiar states u_„, note that the peculiar ground state u_i with n ~ 1 has eigenenergy (mass) 



i\J2 + 2/^1 + a2/(l/2 - ^/l/^^^Y - V2m^/lTa 



«;_! = m Y 2 + 2/\J 1 + ^2/(1/2 - Vl/4 - a^? ^ V2mVl + a, (39) 

which represents very tight binding, with a binding energy on the order 300 keV for an e"*'e~ state and a root-mean- 
square radius on the order of a Compton wave length instead of an angstrom. In particular we find (see Appendix 
C) 



v/pTT^-. (40) 



m 

We note further the anti- intuitive behavior of the peculiar ground state energy (mass), increasing with increasing 
coupling constant a instead of decreasing. The excited states are quite near to the usual bound states. Wc find the 
following pattern for those excited peculiar states 

w-n = 2m - raa^lMji - \f + TOa^/2(n - 1)^ (1 + ll/32(n - 1)) + O(a^); n = 2, 3, 4, ... (41) 

In the nonrelativistic limit, where terms of order a^ are ignored we find that the states are degenerate with the nth 
usual state identical to the (n -}- l)th peculiar state. If we include the ol^ corrections then we find that 

w+„ - W-(„+i) = ~ma^ jr?. (42) 

For all of the usual states and the remaining peculiar states we have 

(r2)+„ ^ i^±^±l^[(n + A+)2 + W + 3], n = 1, 2, 3..., 
6(£u,+„a) 

2\ - ^" + ^"^2[(n-HA_)2 + 5a2 + 3],n = 2,3..., (43) 



r 



6 (£«,_„")" 



so that the size in the {n + l)th peculiar state is nearly the same as with the nth usual state. 

As shown in the Appendix C, the two sets of solutions, are not orthogonal with respect to one another. For 
example, the two n = 1 wave functions have the respective forms 

U-i_(r) = c^r ++"'" exp(— K_)_e^,^ar), 
2 1 



1 + VI - 4a2 A+ + 1 
U-{r) ^ c^r "^"'" exp(— K_eu,_ar), 

2 1 

K- = , = -. , (44 

1 - Vl - 4a2 A_-^l' ^ ' 
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where for brevity of notation, we have omitted the principal quantum number designation in u± for the of the ground 
state. Clearly since they are both zero node solutions we have 



{u-\u+) 



dru+(r)u-(r) ^ 0. 



(45) 



How do we reconcile this with the expected orthogonality of the eigenfunctions of a self-adjoint operator corresponding 
to different eigenvalues. In the present context, the naive self-adjoint property requires that 



(«+l(- 



a 



dx^ 



)W-) = {u-\ 



dx'^ 



-)\u, 



This boils down to 



dxuA 



(Pu_ 
dx"^ 



dxu^ 



Let us integrate by parts. Then we have 



dxuA 



dx-^ 



u+ 



du- 
dx 

du-\- 
dx 



d-'u^ 
dx-^ 

du+ du- 



dx- 



dx dx 
du-L- du^ 



(46) 



(47) 



(48) 



Q dx dx 



We thus have a self-adjoint operator if 



u+- 



du- 
dx 



du^ 
dx 



(49) 



Now clearly these vanish at the upper end points. Since we have that 

1 



du+ A+ + 1 

da; X 



A+ + 1 
1 



du_ _ A_ + 1 
dx X A_ + 1 



(50) 



at the lower end point the LHS of Eq. (|49|) is 

A_ + l 



lima;^++^ exp(-a;/(A+ + l))x 

a:— fO 



A_+l 



exp(-a:;/(A- + 1))( 



A_ + l 



A_ + l- 



limx( 



A_ + l 



(51) 



whereas at the lower end point the RHS is A-|_ + 1 7^ A_ + 1. Thus, the second derivative is not self-adjoint in this 
context! This accounts for the non-orthogonality of the usual and peculiar ground states in Eq. psj) . 
In general beginning with a set of usual and peculiar wave functions {w+n,w_„} such that 



''+n|"+n' 



\^ — n|W — n'/ ^nn' ^ 



(52) 



we find with 



H = 



h±r. 



where x = Su,ar, that 



1 


\ "^^ 


2ewa c? 




\ '^^ 


2 


(swaf 


dr^ 


r r2_ 




dx^ 


X 


h±nU±n, 








-l^±n = -l/(^± 


+ nf, n = 1,2, 






{u+n\H\u+n') = 


0n7i' i^+n 






(u_„|iJ|u_„') =: 


Ojijif ll — ji^ 






{U-n\H 


\u+„') == 


Ojiji' il-\-.fi' 









(53) 



{u+n\H\U-„') = bnn'h-n- ^ (u_„' |iJ|u+„) . 



(54) 
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In the first two terms it docs not matter whether the H operators operate to the left or the right. In the last two 
cases we explicitly have H operating to the right. To emphasize that we write them as 



{U-n\{H\u+n')) = hnn'h+n', 
(u+„|(i/|u_„')) = bnn'h-n'- 



(55) 



It is evident that with both sets of basis, H is not self-adjoint since (w_„|(iJ|u+„')) ^ ((ii_„|iJ)|'u+„) and 

{u+n\{H\U-n')) ^ [{u+n\'H-)\U-n') ■ 

Let us see where the non-orthogonality leads us if we treat both basis on an equal footing. In that case a general 
wave function for the ^5*0 system would be expanded a^ 



* = ^ c+„u+„ + ^ c_. 



n+ 



and applying the variational principle to 



(i/) = 






(56) 



(57) 



and defining (we show just a finite n x n portion of the matrices) 



B 



H. 



H 



bii bi2 . 

bnl bn2 ■ 

h+i 
h+2 

. 
7i_i 

/l_2 




then in block form we would have the eigenvalues equation 

H+ BH 
BH, H 



bin 
b2n 

Onn. 

■ 



h+n. 

■ 




(58) 





c 


= -k' 


1 B 
B 1 




c 



(59) 



(Note that the way this stands , the matrix on the left is not self-adjoint.) Multiplying both sides on the right by 



1 B 
B 1 



(l-B^)-i -B(l-B^ 



-B(l-B^) 



2a-i 



B^) 



2^-l 



we obtain 



H+ 



c 



c 



H^ ^ ^ 

It is clear that the eigenvectors corresponding to the eigenvalue sets of — k1 and — k^ are of the form 








c 



(60) 



(61) 



(62) 



From Eq. (|36p . one recalls that the two sets of basis functions {u+„,m_„} have distinctly different behaviors at the 
origin, corresponding to the usual and peculiar solutions. In particular 

U+nix) = C+nX^^ + ^eM~x)Lll+ + ^{x), 

u_n{x) = c+nx^-+^ exp{-x)Lll-+\x), (63) 



® Strictly speaking we should include the continuum states. See section below for discussion of those states. For the purpose here the 
use of discrete states is sufficient. 
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These generalized Laguerre polynomials are orthonormal with respect to different weight functions a;^±"^^ exp(— a;). 
They would each correspond to a complete set. Together they would constitute an over-complete set. However, that 
does not imply that Eq. (|56p is incorrect as it allows for the function ^(x) to be a linear combination of functions 
with two distinct behaviors at the origin. Nevertheless, the set up here is a bit clumsy with questions of completeness 
and the non-self-adjoint property remaining. 

It should be realized that for the given quasipotential of the type —a^/r"^ at short distances that is at hand, both 
the set of usual states and the peculiar states are physically admissible states. There does not appear to be reasons to 
exclude one set as being unphysical, if one is given the attractive interaction near the origin as it is. We note however 
that the peculiar states with the r~^ behavior at the origin are excluded from existence if coefficient A(A + 1) for the 
l/r^ term is greater than zero since that would lead to a u{r) that is singular at the origin. Only for interactions with 
sufficient attraction at the origin (so that —1/4 < A(A + 1) < 0) can these states be pulled into existence and appear 
as eigenstates in the physically acceptable sheet, with regular non-singular radial wave functions at the origin. It is 
desirable to find ways to admit both types of physical states into a larger Hilbert space to accommodate both sets of 
states with the mass operator to be self-adjoint and the states to be part of a complete set. It is reasonable to assign 
a quantum number which we call "peculiarity" for a states emerging into the physical sheet in this way as physically 
acceptable states. The introduction of the peculiarity quantum number enlarges the Hilbert space, allows the mass 
operator to be self-adjoint, and the set of physically allowed states become a complete set, as we shall demonstrate. 

We introduce a new peculiarity observable C with the quantum number peculiarity <^ such that 

C,X+ = C,X+ with eigenvalue C = +1, 

Cx- = CX- with eigenvalue (^ = — 1, (64) 

with the corresponding spinor wave function xc assigned to the states so that a usual state is represented by the 
peculiarity spinor x+i 

X+ - Q , (65) 

and a peculiar state is represented by the peculiarity spinor x-, 

X- = (?) . (66) 

With this introduction, a general wave function can be expanded in terms of the complete set of basis functions 
{u_i_„, u_„} as 

'^ = ^ainUc,nXi, (67) 

where n represent all the spin and spatial quantum numbers of the state and Q the peculiarity quantum number. The 
variational principle applied to 



would lead to 



Hu+nX+ = -'«+„'"+nXH 



Hu^nX^ ^-K,_^U^nX~- (69) 

It is clear that in this context the usual and peculiar wave functions are orthogonal, H is self-adjoint, and the basis 
states are complete. That is, 



/"OO 

Jo 



(70) 



and so the set of basis functions {ii+„C+7 ""-nC-l; containing both the usual states and peculiar states in the enlarged 
Hilbert space, form a complete set. We also have 



{i\H\j) = {uQ^n^H\uQyn,) = / druc^^n^XQ^HuQ^n.XCj 

Jo 

= (jWN), (71) 
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so that the mass operator H in this enlarged Hilbert space is self-adjoint. We see that the introduction of the 
peculiarity quantum number resolves the problem of over-completeness property of the basis states and the non-self- 
adjoint property of the mass operator. 

2. So Scattering States 
The ^5*0 state equation 

'—^ ^>u = bu (72) 



' 2™a,L(L + l),_2^^^^^^^,^^ ^^3^ 



has the same form as the nonrclativistic Schrodinger equation for Coulomb interaction 

d^ 

except that the standard angular momentum term with L{L + 1) now take on the value of (— a^). The two solutions 
the above equation are given by the regular Fl and irregular Gl Coulomb wave functions, 

u ~ aFL{r],kr) + cGl{ii, kr), 

^7 = - — , (74) 

with only the regular Coulomb wave function having an acceptable behavior at the origin. The long distance behaviors 
of the regular and irregular solutions are 

Fi(?7, fcr — > 00) — ?> const x sin(fcr — rilog2kr + (Tl — Ltt/2), 

GLiv, kr -> 00) -^ const x cos(A:r — rilog2kr + ul — Lit/2), (75) 

in which aL is the Coulomb phase shift given by 

crL = arg(r(L + 1 + ^77). (76) 

Now we can solve Eq. ([72]) exactly for 6^ > by analytically continuing the above solutions to an arbitrary 
(non-integer) angular momentum A and making a few obvious replacements by analogy, 

u = aFx{ri, br) + cGx{r], br), 
A(A + l) = -a2, 

V = ^- (77) 

Using the expressions for the analytically continued Coulomb wave functions to non-integer A [4l| we will presently see 
that we have solutions given by the F and G functions in Eqs. (|83|) and (|84|) below. We emphasize that both solutions 
have an acceptable behavior at the origin. Since A is not an integer, one can replace the irregular solution G\{ri, br) by 
_F'_A-i(?7, 6?')cj In particular, as shown in Appendix D, in terms of the confluent hypergeometric function AI{a, b; z) 

Fxip) = Ca(77)/+i exp(-zp)Af (A + 1 - i?7, 2A + 2; 2ip), (78) 

one has with 

a;(A, 77) = (A + -)tt + a-A-i(??) - <7xiv). (79) 



that 



G,{p) = F-x-i{p)-cos^i\v)FAp) ^ (80) 

sm x{X,ri) 



^^ The reason that Gl is used in place of _F_/,_i for L integer is that the latter is not linearly independent of Fl in that case. It is melded 
together with F^ to produce Gl by a limited process analogous to how the Neumann function is obtained from the Bessel functions. 
For A 7^ integer, F\ and F—x-i are linearly independent. 
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a linear combination of Fx{p) and F_\_i{p). In others words, Eq. (j77|) can be written as 

u = dFx{v,br)+eF^x-i{ri,br), (81) 

where 

A=i(-l + v/l-4a2) = A+, (82) 

-A - 1 = i(-l - Vl-4a2) = A_, 

corresponding to the separate C = il sectors. As with the solutions in Eq. ([TT]). F\{ri,br) and F-\-i{T],br) have 
acceptable behaviors at the origin corresponding to Eq. (|34p. Their respective long distance behaviors are given by 



^a(??, br — > oo) — > const x sin(&r — 7^ log 26?' + a\^ — A+7r/2), 
F-\-i{7], br -^ od) -^ const x sin(6r — r^ log 26?' + ax_ — A_7r/2). (83) 

Alternatively we can use the related G functions to determine the behaviors 

Gx{i],br -> oo) — > const x cos(fer — 77log26r + (Ta+ — -^+^72), 
G-\-i{7],br -^ cx)) -^ const x cos(&r — ?7log26r + ax_ — A_7r/2). (84) 



The respective total Coulomb phase shifts for Eq. ((72|) are the phase shifts for the usual and peculiar solutions over 
and above those due to any angular barrier part (absent here). They are given by 

^A± = cTAi - A±7r/2, 

aA± =arg(r(A± + l + i7?), (85) 



in which 



argr(A± + l + ij])^ r]7P{X± + 1) + ^Z ( x '\ , arctan( 



V 



X± + l + n 'A± + l + r 

(86) 



with the digamma function given by 



^(A± + l) = -7 + A±C(2)-A^^^ (87) 

^ n^(n + A±) 

n— 1 ^ ^ 

The (modified) Coulomb phase shift a\^—X±TT/2 is that for the Coulomb 2ewA plus —A^ term alone. (Again, the ± 
sign corresponds to the two sectors ( = ±1, with usual (+) and peculiar (— ) boundary conditions given in Eq. ([M)) .) 
Without the —A^ term the phase shift would be simply o-q. 



V. SOLUTIONS OF THE TWO BODY DIRAC EQUATIONS FOR THE ^Pq STATE 

A. The ^Po quasipotential 

We now consider the case of the ^Pq state of a fermion-antifermion pair with electric or color charges interacting 
through an electromagnetic- type interaction arising from the exchange of a single photon or gluon. As with the ^Sq 
state, the single photon annihilation diagram does not contribute because the ^Pq state is a charge parity even state. 
Then, the two terms in Eq. ([25]) that precede the V'^A term precisely cancel the barrier term 2/r^ at very short 
distances to give the equation for the radial wave function 

::^b^U. 

The cancellation of terms takes place in the following way. In Eq. (PS|) . the three terms beyond —A^ arise from a 
combination of spin-orbit, spin-spin, tensor and spin-orbit tensor interactions. From a detailed examination of Eq. 



d^ 2 


2ewa 
r 


a^ ^ 87rar(5(r) 
^2 yj^ ^ 2a 


dr2 (r + 2a/w)2 
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(jB12p in the Appendix B, we can see that the spin-orbit and tensor terms gives rise to the first "magnetie interaction" 
term on the right hand side of Eq. (|B12p that has a strongly attractive —^a/wr^ attractive part down to distances on 
the order of 2a /w after which this magnetic interaction approaches — 4/r^. The dominance of the attractive magnetic 
interaction at short distances that can overwhelm the centrifugal barrier is in agreement with the simple intuitive 
classical picture presented in the Introduction. The second term on the right-hand side of Eq. (jB12| . arising from 
a combination of Darwin, spin-spin and tensor terms, has a stronger repulsive %a^ jixP'r'^ part down to distances on 
the order of 2a /w after which it approaches -|-2/r^. Together they tend to exactly cancel the angular momentum 
barrier term +2/r^ at very short distances. In addition to the repulsive interaction containing 6{r) arising from the 
assumption that the electron and positron arc point particles, the quasipotential behaves as —a^/r^ at short-distances, 
separated from the outside long-distance region by a barrier. The interaction containing the delta function comes 
from a combination of Darwin, spin-spin, and tensor terms. Three fourths of the repulsive term containing 5{r) 
comes from the Darwin piece while one fourth from the combination of the spin-spin and tensor parts. For brevity of 
nomenclature we shall just call it the delta function term. 

One of us (HWC) examined in a previous work [40| the effects on bound state energies due to a repulsive 5{r) 
interaction by itself, without additional radial dependence. It was found that for wave functions ip that do not vanish 
at the origin and for potentials that are less singular than 1/r^, the exact effects on the eigenvalue of including a 
repulsive delta function do not agree with the results of perturbation theory in the limit of weak coupling, when the 
delta function potential is modeled as the limit of a sequence of spherically symmetric square wells. In particular it 
is shown that the repulsive delta function, viewed as the limit of square well potentials, produces no effects at all on 
bound state energies. In our case here the appearance of the 5{r) potential differs from this reference in two aspects 
however. First of all the 5{r) appears in conjunction with r/{'wr + 2a), softening its repulsive effects. Secondly, the 
wave function ip = u/r for the solution without the delta function term diverges at the origin both for what we call 
the usual solution and what we call the peculiar solution. If the null effects on bound state energies and phase shifts 
seen in [40| should occur in our case as well, this, however, does not lead to a problem with perturbative agreement 
with the spectral results. 

In the case of weak potentials where the denominator (wr + 2a) is replaced by wr, we have shown previously in 
[la | that the remaining terms in Eq. (j88|) without the delta function term, when treated nonperturbatively, would 
produce numerically the same spectral results for the '^Pq state as the inclusion of the repulsive S{r) interaction treated 
perturbatively. The agreement of the perturbative treatment with the delta function term for weak coupling with the 
nonperturbative treatment containing no delta function term justifies the first approximate analysis of ignoring the 
delta function term and treating the remainder of the equation nonperturbatively in the following subsection. 

B. Usual and Peculiar Solutions for the "^Pq State 
The wave equation ([55|) for the ^Pq state without the delta function term becomes 



(P 2 2e,„a -^ 



a I ,2 



dr"^ (r + 2a/w)'^ r r^ 

with a short distance (r << 2a /w) form 



-Tfu^ b'u, (89) 






-— --^ = 0, (90) 



the same as with the ^5*0 states. Thus, the ^Pq states also have the same types of solutions as the ^Sq states, with 
radial wave functions near the origin as given in Eqs. ([33]) -([34]). Thus, there are usual '^Pq states with peculiarity 1, 
and peculiar ^Pq states with peculiarity —1. 

Note that both the usual and the peculiar solutions u± ^ j,(i±\/i-4a^)/2 ^j-jg^ fj-om the strong magnetic interaction 
that significantly modifies the qualitative behavior of the interaction at short distances, when the total spin and the 
orbital angular momentum are oppositely aligned in the ^Pq state. If the strong magnetic interaction is absent, the 
2/(r + 2a/w)'^ term in Eq. (|55)) would be 2/r^, and the wave function near the origin would be 

with 

V-^d-V = r^<^'^^^^^^^^^drdn. (92) 

In that case, as stated below Eq. ([35)1 . only the usual u+ solution is quantum- mechanically admissible, while the U- 
state becomes singular at short distances. Such a comparison shows that the peculiar solution u„ is not present when 
there is no strongly attractive magnetic interaction at short distances or more generally for J ^ 0. 
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C. The 5 function term and the charge distribution 

The discussions in the above subsection pertain to the quasipotential without the delta function term. We now 
examine the full Eq. (|25l) for both the usual and peculiar solutions with the 5{r) term included. Consider first the 
perturbative treatment of taking the interaction containing S{r) as a perturbation. We evaluate the expectation value 
of the interaction term containing S{r). Even though both usual and peculiar solutions have a diverging ^±{r) near 
the origin they each are allowed as a probability amplitude since the probability J^y |'0±(t)| cfir for an arbitrarily 
small volume AV about the origin would be finite, in addition to the essential boundary condition u±(0) = 0. With 
|^±(r)| near the origin having the behavior of r'^^^^^'^^^\ the expectation value oi d{r)/{w — 2A), after performing 
the angular integration, is 



J wr + 2a J 



wr + 2a 



5{r) 



{wr + 2a) ^ ' 

which is zero for the plus sign for the usual solution but diverges for the minus sign for the peculiar solution. 

The results of Eq. (j93p for the usual solution explains our previous agreement between (i) the perturbative treatment 
with the delta function term for weak coupling and (ii) the nonperturbative treatment without the delta function 
term |16l | . The agreement arises because in Ref. [l6[ we limited our attention only to the usual solution for which the 
expectation value of the delta function term is zero. 

The results of Eq. ([M|) for the peculiar solution indicates that the delta function term cannot be treated as a 
perturbation in the present formulation, as such a treatment will lead to a diverging energy. The delta function term 
arises from the charge distribution of the interacting particles, as it is related to the Laplacian of the gauge field, V^ A, 
as given in Eq. ([26]). A proper non-perturbative treatment of the problem of the peculiar solution states requires 
the knowledge of the wave function at very short distances. Therefore, it will require not only the knowledge of the 
structure of the charge distribution but also the necessary auxiliary interactions at even shorter distances that are 
needed to bind the charge elements of the distribution together. The auxiliary interactions will affect the solutions of 
the two-body wave functions at very short distances and the states of the peculiar solution. At the present moment, 
we have little knowledge of the structure of elementary charges, much less the auxiliary forces that would bind the 
charge distribution together at very short distances. 

The structure of the charge distribution of elementary particles at very short distances is basically an experimental 
question. As the strong magnetic interaction allows the two interacting particles to probe the short-distance region, 
it is therefore useful to investigate quantities that may reveal information on the structure of the charge distribution. 
While many possibilities can be opened for examination, we shall examine the following possibilities in the present 
manuscript: 

(i) We shall first examine the case in which the (unknown) auxiliary interaction that binds the charge elements 
of the elementary particles together and the repulsive interaction arising from the charge density p{r) counteract in 
such a way that the total interaction at short distances would still be dominated by the —a^jr^ term. Under such a 
circumstance, the effects of the auxiliary interaction would cause the delta function term term in Eq. (j88p to make no 
contribution at short distances. Keeping the dominant terms, the equation of motion for the wave function becomes 
Eq. ([55]) without the delta function term. It also must be recognized that for the usual solution, the perturbative efi'ect 
of the delta function term (in which we ignore the effect or the potential in the denominator w — 2A) is accounted for 
by a nonperturbative (numerical) treatment of the entire $ without the delta function term. So, our treatment of 
the delta function term in this case parallels that used in our earlier spectroscopy calculations [l6| . 

(ii) We examine subsequently the case when the auxiliary interaction that holds the charge element together leaves 
the gauge field A(r) unchanged while the delta function term in Eq. (|88p is modified by treating the delta function 
as the limit of a set of Gaussian distributions with different widths. 

(iii) We examine two additional models completely within QED (or QCD) with an assumed basic charge distribution 
that generates the gauge field also in the region interior to the charge distribution. However, the auxiliary interactions 
that hold the charge together and that can interact with the other antifermion are altogether neglected. It should 
be recognized that within pure QED (or QCD), with the neglect of the auxiliary interactions that hold the charge 
elements together, the charge distribution cannot be a stable configuration. 

In the next section we describe the method we use to indicate the presence or absence of a resonance in the "^Pq 
system. 
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VI. PHASE SHIFT ANALYSIS 

111 our study of the "^Pq state for both the usual and pecuhar solutions, we wish to find out whether or not there is 
an energy that will lead to a 7r/2 phase shift for a given a and constituent mass m . Equation (P5|) for the "^Pq state 
is a Schrodinger-like equation of the form 

-^ + ^^^^^ ^^ + ^{r)\u{r) ^ bMr)- (94) 

We calculate the phase shift for this problem by the variable phase method of Calogero [l^ . We first describe this 
method generally (see Appendix E for a more detailed review) and then later in this section describe its application 
to the '^Pq state. We take W{r) to include not only the quasipotential $(r) but also the angular momentum barrier. 

W{r) = ^^^^^ ^^ + $(r). (95) 

Thus our equation has the form 

j^ + W{r) ^ u{r) ^ b^u{r). (96) 

The Calogero method relies on introducing a reference potential W{r) that can be solved exactly, with two independent 
solutions ui and U2, 

(P - ^ 
-^ + W{r) \ u,{r) = b^u,{r), z = 1, 2. (97) 

There are many ways to choose the reference potential W{r). To display the general idea, we consider the case in 
which W{r) is short range. In that case the phase shift 6l is defined by 

u{r-^oo)^sm{br-LTr/2 + 6L)- (98) 

The Calogero method uses two different types of W{r). In the first, W{r) = Wj{r), the reference potential has the 
same long and short distance behavior as W{r). In the second W{r) = Wii{r), the reference potential does not have 
the same long and short distance behavior as W{r) but is especially simple. 

We consider first Type I reference potential, Wi{r) = L(L + l)/r^, the angular momentum barrier potential, for 
which the reference wave functions Mi(r) and U2{r) are the well known spherical Bessel functions jL(&r) and n]^{br) 
[l9|. The solution ui{r) is taken to be the regular solution, having the same short distance behavior as u{r), in 
particular, ■ui{r — >■ 0) = 0. The solution U2{r) is taken to be the irregular solution, U2{r — >■ 0) ^ 0. Those functions 
together with their long distance behaviors are given by 

ui(r) = jL{br) — > const sin(fer — Lt:/2), 

U2{''') = —nL{br) — > const cos(fer — Ltt/2). (99) 

We introduce the amplitude function a(r) and phase shift function (5L(r) to represent the wave function solutions 
for the W{r) potential, u{r) and u'{r), as 

u{r) = Q!(r)(cos SL{r)ui{r) + sin(5i(r)-U2(r)), 
u'{r) = a{r){cosSL{r)u[{r) + sin SL{r)u2{r)). (100) 

This leads to the following equation for the phase shift function (see Appendix E) 

[U2(r)u'{r) — U2{r)u(r)) 
Further manipulations lead to the differential equation for 6L{r) given by 



S'Lir) = - 



[Wjr) - W{r)] [. w. A • A ^ aI' 
jL{br)cosdL{r)-nL{br)smdL[r) , 



(102) 
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To find the connection to the phase shift (5^ note that from Eq. ( I100[) and ([M|) . 

u(r — ?► oo) = const {cos (5^ (r ^ oo) sin(6r — Lii/2) + sin 5l {r — > oo) cos(6r — LTr/2)} 

= const X sin(&7 — L'k/2 + 5l{oo)), (103) 

and so comparison with ((M|) gives the solution of the phase shift 6l for the W{r) potential as 

6l = Sl{(x). (104) 

Thus, the second order linear different equation becomes a first order non-linear equation whose solution at ?' — > oo 
gives the phase shifts of the scattering problem with the W{r) effective potential. The boundary condition of (5^(0) = 
follows from Eq. (|10ip when one chooses ui{r) to have the same behavior as u{r) as r — ?> 0. 

We consider next type II of the short-range reference potentials Wu (r) which do not need to have the same long 
distance behavior as W{r) as long as the Schrodinger Eq. ((97|) containing the reference potential Wii{r) has an exact 
solution. For example we may choose Wii{r) = 0. Then the two exact reference solutions of Eq. (P7| are simply 

ui{r) = sin(6r), 

U2{r) = cos(6r). (105) 

One defines a phase shift function jLir) as in equation ()100|) so that 

u{r — > oo) = const x {cos7i(r — > oo) sin(6r) + sin7i(r — >■ oo) cos(6r)} 

= const X sin(6r + 7i(oo)). (106) 

Comparison with ((98)) gives 

'5l-7l(oo) + ^- (107) 

Since the angular momentum barrier is excluded from the equations for Ui (r) one finds that the phase shift equation 
for integrating the phase shift function 7l('^) includes the repulsive barrier term in W [Eq. ([M]) ]. 






b 



W{r) ^., 



cos7L(r) sin(6r) -I- sin 7^(7-) cos(6r) 
sin2(fer + 7i(r)). (108) 



b 

Note that because of the L{L + l)/r^ behavior of W{r) — Wii{r)(r) = W{r), which dominates at large distances, one 
will have to integrate quite far to obtain convergence for 7L(r)13 For this case of Wii{r) = 0, one has an equation 
similar to (|10ip with SL{r) replaced by 7L(r). Thus even though ui{r) has a different behavior than u{r), we still 
have the boundary condition 7^(0) = 0. Eq. (|107p compensates for the —LTr/2 effective phase shift due to the barrier 
term in W{r) in Eq. ( llOSp . We also have the additional boundary condition of (see Appendix E) 

7i(0) = -^. (109) 

We now turn our attention to the ^Pq system, in particular Eq. ([^5]) for a general $(r). In this application of 
the Calogero method we choose a reference potential W{r) = Wiii{r) that in a sense is a hybrid of the two types of 
reference potentials considered above. Since Eq. p8p contains a long range Coulomb interaction — 2£u,a/r we must 
include that interaction into our choice for Wiii{r). If it did not have the same behavior as W{r) at large distances 
we would have to have a way of subtracting an infinite Coulomb phase shift, log 2br . So, in this way our application 
is similar to the first type Wj (r)above. We also include the ~c? jr^ term in Wiii(t) because as seen in Eqs. (|89p . 
(|M| and ([M)) the solution displays the desired short distance peculiar as well as usual behaviors. We do not include 
the angular momentum barrier term 2/r^ however, as this would prohibit a treatment of the peculiar solution (see 
comments below Eq. (|M|)). Thus we choose 

WuAr) = -, 

W{r) = \^^{r), (110) 



Alternatively Calogero gives a formula for avoiding integrating to large distances to build up a centrifugal phase shift. (See [l9|, p 92). 
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where $(r) is given in Eq. ([2^. In this way Wiii{r) has some similarities to the second type Wij{r) discussed above. 
Our choice for W7//(r) permits the two exact solutions ui{r),U2{r) of Eq. (|97p which becomes that of the ^5*0 state 

■^ 7^ -^\ - = ''-■ (111) 

Now to determine the phase shift for the actual '^Pq state we return to the conditions defined in Eq. (|110p . Then the 
full solution has the asymptotic form 

u{r — )• cx)) -^ const x sin(6r — ri\og2br + cti — tt/2 + 5i). (112) 

The appearance of ai and i5i includes the effects of the angular momentum barrier term 2/r^ in the presence of the 
Coulomb interaction. In Appendix E, using 

wi(r) :=: Fx{r],br), 

U2{r)^Gxiv,br), (113) 

we show that the full ^Pq phase shift 5 is given by 

6 = 5, + a,= 7±(oo) + aA± + (1 - A±)^/2, (114) 

where (in analogy to the proof of Eq. (|108p with W ^ 0) 7±(r)satisfies the nonlinear equation 



l±ir) = 



W{r)-Wiiiir) \ . . . , ,^ '^ 
^^-— ^ cos 7± (r)i^A± +sm7±(r)GA± 



(115) 



subject to the boundary condition that 7±(0) = (see Appendix E). The functions F\j_ and G\^ are the regular and 
irregular Coulomb wave functions corresponding to the negative effective centrifugal barrier —ct^jr^. Again, because 
of the 2/r^ behavior of W(r) — Wmir) which takes over at large distances, one will have to integrate quite far to 
obtain convergence for 7± (r) . 

We consider numerical solutions for both the usual solution with A+ = (— 1 + Vl — 4a2)/2, and the peculiar solution, 
with A- = ( — 1 — \/l — 4a^)/2. In the next section we discuss the results obtained in the numerical integration of the 
phase shift equation ( I115P for different behaviors at very short distances. 

VII. NUMERICAL RESULTS FOR ^Po RESONANCES 

A. The case without the delta function term 

With the above general formalism, we can begin to examine states in the quasipotential of Eq. (|88| first without 
the delta function term. The Schrodinger equation for the ^Pq state becomes Eq. ([M|) . In order to gain an idea on the 
attractive magnetic interaction at short distances for this '^Pq state, we plot in Fig. 2 the corresponding quasipotential 
including the angular momentum barrier, 

W{t) = ? 26^ _ a^^ 

(r + 2a/ w)'^ r r^ ' 

for w = 27.85 MeV, a ~ 1/137 and a constituent electron mass of 0.511 MeV. One observes that at short distances 
W{r) becomes very attractive and behaves as —a^/r^. There is a barrier in the region between 10^^ to 10^^ GeV^^. 
Such a potential becomes singular at r — > when a exceeds 1/2 [1|. 

We calculate the phase shift as a function of energy using the boundary condition 7±(0)=0, including the dependence 
of the potential as a function of energy. For the usual solution [C, = +1), our results for the QED e~e^ system in 
the ■^Po state with a — 1/137 and m ~ 0.511 MeV show no evidence whatsoever for resonances for all cm. energies 
tested (from about 1 MeV to about 100 MeV). The magnitude of the phase shifts are of the order of tt/IOO. 

For the peculiar solution (C — —1) with the wave functions starting with a less positive slope, the attraction at 
short distances is able to bend the wave function downward to result in a very sharp resonance at about 27.85 MeV. 
In Figure 3(a) we plot the phase shift 5 ~ 5i+ai as a function of the cm. energy w and sin^ 5 versus w in Fig. 3 (b). 
We start the integration at the origin and extend to about 1 angstrom. As one observes, the phase shift undergoes a 
transition from near zero to tt. The resonance has a full width at half maximum of 15 KeV. We also include a plot of 
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FIG. 2. The effective potential W{r) = 2/{r + 2a/wf' - la^ajr - d^ jr^ for the for the (e+e") system in the ^Po state with 
a = 1/137 and w = 27.85 MeV. 

TABLE I. Variation of the resonant energy as a function of the quark mass for a fixed ol^ = 0.11. 



quark 


mass 


■Wr 


up 


3 MeV 


27 MeV 


down 


5 MeV 


45 MeV 


strange 


135 MeV 


1220 MeV 


charm 


1.5 GeV 


13.6 GeV 


bottom 


4.5 GeV 


40.8 GeV 


top 


175 GeV 


1590 GeV 



the wave function in Fig. 4 from the origin up to about 1000 GeV~^. The wave function rises as j-(i-vi-4a )/2 j^ga,r 
the origin, and appears nearly flat at r ^ 10~^ GeV~^, and it slowly decreases near the barrier. It oscillates when it 
emerges from the barrier atr^2xl0~^ GcV~^. 

Having observed a resonance for the QED interaction with the e+ and e~ constituents, we turn our attention to 
quarks and antiquarks interacting with a color-coulomb type interaction with an effective coupling constant a^. We 
focus here only on the ^ = — 1 sector. In the color-singlet {qq) states of interest, the effective interaction is then 
Ooff = 4as/3. To get an idea of the order of energy for these quark-antiquark two-body resonance states, we calculate 
the resonance energies for the typical case of ots = 0.11 For this value, the resonance energy varies nearly linearly with 
quark mass. The largest energy resonances occur with the largest quark masses. In Table I we present the resonance 
energies w^ for the families of quarks from the up quark to the top quark. It should be pointed out that these 
resonance values take into account only the Coulomb-like portion of A{r) = — (4/3)as/r and ignores any affects on 
the resonance values of the confining part of the potential. 

To examine how the resonance energies varies with the coupling constant, we have found that for fixed mass (e.g. 
0.511 MeV) the resonance energy wr increases as the coupling parameter decreases until the coupling constant gets 
to be on the order of 0.01, when wr starts decreasing again. 



B. The case of representing the delta function by a Gaussian function 

For the second case for the ^Pq state given in Eqs. ([55]) and f 1^5)) using (jllSp . we take 



W{r) = 



2e„a a^ 87rar(5(r) 



(r -|- lajw)"^ r r^ [wr + 2a) ' 

in which we model the three dimensional delta function by 

exp(-rV2cr^) 



(5(r) -> (5^(r) = 



(271)3/2^3 



(116) 



(117) 



22 



4 




1 


1 




1 


3 

«>2 

1 






(a) 








- 






: 


- 


^^ 






1 












1 

0.8 
:0.6 


- 


(b) 




(e 
a 
w 


' 1 

e ) system — 

= 1/137 

= 27.85 MeV - 


0.4 


- 








- 


0.2 


- 


^1 






", ^ : 


82 


27.84 


27.86 
w (MeV) 




27.88 



FIG. 3. The phase shift as a function of w for the (e^e ) system with a = 1/137 and w = 27.85 MeV. 
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FIG. 4. The wave function it of the peculiar resonance at to = 27.85 MeV for a — 1/137 and m = 0.511 MeV. 



In this treatment, we keep the point charge source term for the A{r) so that A{r) = —a/r. What we are attempting 
to do is just present a mathematical representation of the delta function that will allow a numerical solution. The 
reference potential Wiii{r) is the same as without the delta function. With this modeling of the delta function we 
start off our Rungc-Kutta integration of Eq. (|115p with 7-(0) = since W{r) — Wiii{r) is w'^2a^ at the origin just 
as without the delta function term. The function 6^ docs not alter the extreme short distance behavior since it is 
multiplied by r and vanishes at the origin. Wc obtain the resonance energy results as given in Table IT It is obvious 
that for small tq wc obtain a limiting behavior oi w = 3.13 GeV-fm/\/2cr. There is however, a difference between 
what we are doing here and what was done in [40|. There the delta function was just regarded as given, not related to 
other parts of the potential. Here that is not the case. The delta function arose from the Laplacian of A{r). There 
may therefore be some ambiguity of, in effect, modeling V^A in one part of the quasipotential while leaving A{r) 
unaffected in the other part. That leads us then to the third case. 
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TABLE II. Variation of resonance energy with the width of the Gaussian distribution 



\/2o-(fm) 


«;(GeV) 


1000 


0.0279 


100 


0.0278 


10 


0.0279 


1 


0.0398 


0.1 


0.314 


0.01 


3.13 


0.001 


31.3 


0.0001 


313 


0.00001 


3130 


0.000001 


31300 


0.0000001 


313000 



C. The case of representing the charge distribution by a continuous function 

In Eq. (|25|) . if one replaces A{r) by 



I a a\ 1 a , . 

^= ]-. n 77^^--' (118) 

\r ro J 1 + cxp|(r - ro)/dro| r 



or alternatively as 



I or- _ 3a_ for r < rn 

Air) = ^^ ^- ' - 'V (119) 

[-f forr>ro. 

then our numerical solutions show that there is no ^Pq resonance for the peculiar solution for both cases, resulting 
in a phase shift of n all the way down to threshold {w = 2m). Both of these corresponds to smeared charge 
distributions from V^ A but neither have auxiliary interactions at short distances that would bind the elements of the 
charge distribution together. The reason no resonance is produced in this case is that in the interior of the charge 
distribution (r < ro), the angular momentum barrier in Eq. ()25|) comes out from under the dominance of the magnetic 
interaction terms as A{r) tends to a finite constant. By following steps similar to those used to determine 7-(0) in 
Appendix E for the point charge one can show this results in an initial value for 7_ (0) defined by 

tan7_(0) = tana;(A,ry). (120) 

This is positive and even though small (^ 0.007) is large enough to prevent the formation of a resonance. Note that 
this differs from the previous section in that here we are giving a physical connection V^A ~ 4Trap{r) between the 
smeared delta function and the invariant potential A{r), whereas in the previous section we simply mathematically 
modeled the delta function in isolation. 

VIII. DISCUSSION AND CONCLUSION 

Magnetic interactions in the ^Sq and ^Pq states are very attractive and singular at short distances. In the two- 
body Dirac equation formulated in constraint dynamics, the magnetic interactions lead to quasipotentials that behave 
as —cP'lr^ near the origin and admit two different types of states. At short distances, the radial wave functions 
u(r) of the usual states, grow as r^+^, while the radial wave functions of the peculiar states grow as r~^, where 
A = (—1 + Vl — 4q;^)/2. They have drastically different properties. 

The existence of usual and peculiar states for the same fermion-antifermion system poses conceptual and mathe- 
matical problems. If we keep both sets of states in the same Hilbert space, then each set is complete by itself, but 
the two sets of states are not orthogonal to each other. Our system is thus over-complete. Furthermore, the matrix 
element of tl (the scaled invariant mass operator for these states) between states of one type and state of the other 
type are not symmetric and the PI operator is not self-adjoint. 

Given our quasipotential of the type —o^jr^ at short distances for the ^S'o and "^Pq states, both the usual and 
peculiar states are physically admissible. There do not appear to be compelling reasons to exclude one of the two sets 
as being unphysical, if one is given the attractive interaction —a'^/r'^ near the origin as it is. It is desirable to find 
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ways to admit both types of states as physical states while maintaining the self-adjoint property of the mass operator 
and the completeness property of the set of basis states. 

We are therefore motivated to introduce a quantum number C, which we call "peculiarity" , to specify the usual 
or peculiar properties of a state. The peculiarity quantum number C is 1 for usual states which have properties the 
same as those one usually encounters in QED and QCD. The peculiarity quantum number is —1 for peculiar states 
which intrudes into the physical region, when the interaction near the origin becomes very attractive, such as the 
A(A + l)/r^ interaction with —1/4 < A(A + 1) < 0. The introduction of the peculiarity quantum number enlarges 
the Hilbert space, makes the mass operator self-adjoint, and the enlarged physical basis states containing both usual 
and peculiar states in a complete set. It is also clear from our discussions that to maintain the self-adjoint property 
of the mass operator and to have a single complete set, the presence of the peculiarity quantum number will be a 
general phenomenon, when the mass operator contains very attractive interactions at short distances such that there 
are more than one set of eigenstates satisfying the boundary conditions at the origin. 

It should be emphasized that the quasipotential —a^/r^ has been obtained under the assumption of a point fermion 
and a point antifermion for which the gauge field potential between them is A{r) = —a/r. The point nature of an 
electron may be a good experimental concept as the lower limits on the QED cut-off parameter Acut with the present 
day high energy accelerators exceeds the value of 250 GeV, suggesting that the electron, muon, and tauon, behave as 
point particles down to 10~^ fm. The asymptotic freedom is a good description for the interaction of quarks at short 
distances. It may appear that point charge particles may be a reasonable description. On the other hand, a finite 
structure of the electron or quarks may modify significantly the short-distance attractive interactions so substantially 
that the peculiar states may be pushed out of existence. The experimental search of the peculiar states, which follows 
from the point charge potential, can provide a probe of the point nature of these particles and the interaction at short 
distances. 

Our first focus on the attractive magnetic interaction is for the ^5*0 states, where the spins of the fermion and 
antifermion of opposite electric or color charges are oppositely aligned. The usual bound ^Sq states possess attributes 
the same as those one usually encounters in QED and QCD, with bound state energies explicitly agreeing with the 
standard perturbative results through order a^. In contrast, the peculiar bound ^5*0 states, yet to be observed, not 
only have different behaviors at the origin, but also distinctly different bound state properties (and scattering phase 
shifts). For the peculiar ^5*0 ground state of a fermion-antifermion pair with fermion rest mass ra, the root- mean- 
square radius is approximately 1/to, binding energies approximately (2— \/2)m^ and a rest mass approximately \f2m. 
On the other hand, the {n -f l)"'^5o peculiar state with principal quantum number (n -I- 1) is nearly degenerate in 
energy and approximately equal in size with the n^S^ usual states. 

Our second focus is for the '^Pq state where the total spin and the orbital angular momentum are oppositely 
aligned. The magnetic interaction overwhelms the centrifugal repulsion at short distances and the wave function 
admits a peculiar solution that grows with radial distances as u ~ ^(i-vi-4q )/2^ "Wvs. particle charge density p(f) 
and auxiliary interactions that bind the charge elements together can be exposed for scrutiny. As the structures 
of elementary particles are basically experimental questions, it is useful to utilize the magnetic interaction to probe 
such charge distributions at very short distances. While many possibilities can be opened for examination, we have 
investigated only a few possibilities in the present manuscript. 

The ^Pq quasipotential contains a term proportional to 8(r) . As the delta function term does not contribute to the 
usual QED ^Pq bound state energies, it was plausible to ignore it as one of our explored possibilities. In that case, we 
find that there is a magnetic ^Pq resonance at 27.85 MeV for the peculiar solution of the (e"''e~) system. For various 
(55) systems of different fiavors, we find magnetic ^Pg resonances at energies of the peculiar solution ranging from 
many tens of MeV to thousands of GeV. It is interesting to note that these ^P^^ resonances have the same quantum 
number as the vacuum. 

In another one of our explored possibilities, if we mathematically model the delta function at short distances by a 
sequence of Gaussians of different widths without changing the gauge field Air) = —a/r, then a completely different 
behavior for the resonance energies ensues as they occur at different energies, depending on the width of the Gaussian. 
In the third of our explored possibilities, if we replace the delta function by a charge distribution that also alters the 
gauge field A{r), we obtain no resonance at all. 

Because of 1) the limited knowledge of the unknown auxiliary interactions and charge distributions at very short 
distances, not to mention possible alterations on the angular momentum barrier itself, and 2) the ambiguity of treating 
the delta function in isolation nonperturbatively, and 3) the fact that the delta function term does not contribute to 
the ^Pq usual bound state solution, we speculate that the first case may provided a more reliable representation of 
the physics. It furthermore makes a clear prediction of a QED resonance in a region that has not been investigated. 

While we have studied the resonance ^Pq states, future work calls for the investigation of possible ^Pq peculiar 
bound states where the attractive interaction near the origin may allow the formation of bound states. The presence 
of a delta function repulsion at the origin will also lead to difficulties and problems similar to the ones we encounter 
here with the '^Pq peculiar resonances. 
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Fermion-antifcrniion states as we know them experimentally belong to the usual states. Peculiar states have not 
been observed. Can the peculiar states be observed? How do the usual and peculiar states interplay between them? 
Will there be transitions between the usual states and peculiar states? Clearly, the stability of peculiar states first 
and foremost depends on the strong attraction near the origin, which in turn depends on the point-like nature of 
the elementary particles. As we discussed earlier, substantial modification of the attractive interaction at the origin 
may push the peculiar states out of existence. Only for interactions with sufficient attraction at the origin can the 
peculiar states be pulled into existence and appear as eigenstates in the physically acceptable sheet, with non-singular 
radial wave functions at the origin. This is true for both ^5*0 and ^Pq states. From such a perspective, we expect 
that interactions at short distances have important bearings on the existence or non-existence of the peculiar states, 
and presumably also on the transition between the usual and peculiar states. However, the interactions at short 
distances that may allow the peculiar states to be stable and may effect transitions between states with different 
peculiarity quantum numbers (flipping the peculiarity spinor) are not yet known. They can only be obtained by 
careful experimental investigations. The first task of such investigations should be to locate these peculiar states in 
high-energy experiments where interactions at short-distance may be involved and these strong interactions at short 
distances may lead us to probe short-distance transition from the usual to the peculiar states. These new ^5'o peculiar 
bound states correspond to a very tightly bound state and a set of (n -I- l)th excited states nearly degenerate with 
the nth usual states. It will also be of interest to search for these states as a result of some tunneling process between 
the usual and peculiar states, relying on the small probability of the usual states to explore short-distance regions 
where the interaction at short distances may induce a transition from a usual state to a peculiar state. The fact that 
peculiar states of {n + l)th^S'o state is nearly degenerate with the usual nth ^Sq state may facilitate such a tunneling 
transition. Whether or not these quantum-mechanically acceptable resonances correspond to physical states remains 
to be further investigated. Future experimental as well as theoretical work on this interesting topic will be of great 
interest in shedding light on the question whether magnetic bound states and resonances play any role in the states 
of fermion-antifermion systems. 

Future work should include the effects of the weak interactions, in particular the exchange of the Z° boson. Since 
the mass of the Z° is about 92.5 GeV the range is on the order of 10~^ GeV~^. The exchange of this particle 
corresponds to not only a vector interaction but also a pseudovector interaction. The coupling corresponding to the 
vector portion is [45j 
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With sin^ 6* - 0.23 we find that 



-0.25e (123) 



so, its coupling appears with the same sign as that of the photon. Since a* = e*^ ^ 0.063 Its effect should be small 
but not negligible. There is also the question of the effects of the pseudovector interaction, not discussed in this 
appendix but in [HI, [H. 

Finally, there are however important mathematical and conceptual issues associated with these two-body fermion- 
antifermion system at short distances that require future careful considerations. In standard QED theory, the charge 
and mass of a single charged object due to vacuum polarization and self energy corrections need to be renormalized or 
regularized to render them finite for comparison with observables. For the case with two-body magnetic bound and 
resonance states, for example, how are the two-body Green's functions regularized, with internal lines off mass shell in 
a way that refiects the Dirac constraints? How do such regularizations modify the short distance two-body interaction? 
Can the regularization affects the magnetic interaction at short distances so substantially that the peculiar states no 
longer survive to intrude into the physical states? Are these peculiar states stable against fluctuation of the vacuum 
in quantum field theory. These are some of the many interesting questions associated with the two-body problem 
raised by the possibility of magnetic states under consideration. 
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Appendix A: Details of the equivalent Relativistic Schrodinger Equation 
1. Connections between TBDE and the equivalent Relativistic Schrodinger equation [Eq. (|17p ] 

Here we present an outline of some details of Eq. (|T4| and its Pauli-Schrodinger reduction given in full elsewhere 
(see |25l. [29l. l4a l43j ) .This appendix and the one following it are specializations of Appendices A and B given in [l4| . 
Each of the two Dirac equations in (|14p has a form similar to a single particle Dirac equation in an external four- vector 
and scalar potential but here acting on sixteen component wave function ^ which is the product of an external part 
being a plane wave eigenstate of P multiplying the internal wave function tp 



^ 



tp2 
V'3 



(Al) 



The four -0^ are each four-component spinor wave functions. To obtain the actual general spin dependent forms 
of those A^ potentials (including scalar interactions in general) which were required by the compatibility condition 
[Si,S2\ip = was a most perplexing problem, involving the discovery of underlying supersymmetries in the case of 
scalar and time-like vector interactions |ll|.|25|. Extending those external potential forms to more general covariant 
interactions necessitated an entirely different approach leading to what is called the hyperbolic form of the TBDE. 
Their most general form for compatible TBDE is 

Si^ = (cosh(A)Si + sinh(A)S'2)V' = 0, 

S2^ = (cosh(A)S'2 + sinh(A)S'i)^ = 0, (A2) 

where A represents any invariant interaction singly or in combination. It has a matrix structure in addition to 
coordinate dependence. Depending on that matrix structure we have either covariant vector, scalar or more general 
covariant tensor interactions |42| . The operators ^i and S2 are auxiliary constraints satisfying 

Sii' = (5io cosh(A) + 520 sinh(A) )V - 0, 

S2tp = {S20 cosh(A) + 5io sinh(A) )^ = 0, (A3) 

in which the Sio are the free Dirac operators 

Sio == ^75i(7i -pi + m,). (A4) 

This, in turn leads to the two compatibility conditions [l^, 120: S^] 

[Si,S2]^ = 0, (A5) 

and 

[^i,^2]V^ = 0, (A6) 

provided that A(a;) = A(a;^). These compatibility conditions do not restrict the gamma matrix structure of A. That 
matrix structure is determined by the type of vertex- vertex structure we wish to incorporate in the interaction. The 
three types of invariant interactions A that was used in the relativistic quark model based on this approach (as most 
recently discussed in [1^ , [ij] ) are 

r f ^ 
^c{x±) = -I1I2 — ^^1' ^1 = -751752, scalar, 

'^j{x±) = /3i/32 — ;:; — Ci, time-like vector, 

^e(a;±) = 7i-L • 72± — 7: — Oi, space-like vector, (A7) 



where 



75i = ihiifii^ 



(A8) 
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For general independent scalar, time-like vector, and space-like vector interactions we have 

A{x^)^Ac + Aj + Ag. (A9) 

The special case of an electromagnetic-like interaction (in the Feynman gauge) applied in this paper and in |16j 
corresponds to J = —G or 

Aj + Ag = Asm = (-71 • P12 ■ P + li±- 12±) — ^ — Oi 

= 71 -72^01. (AlO) 

and for scalar and electromagnetic interaction, 

A{x^) = Ac + A£M- (All) 

This leads tcEl [H, IH 

Si^^{- G(3i^i ■ V2 + Eil3ij5i + Afi75i - G^S2 • 5(/:/32-^/3i)75i752)V' = 0, 
S2^ = (G/32S2 • Vi + E2hl52 + M2752 + G^Ei • a(£/3i-^/32)75i752)V' = 0, (A12) 

in which d^ ~ d/dx^. With 

G = exp5, 

r, = p± - '-j:, ■ dGE,. (A13) 

The connections between what wc call the vertex invariants £, JJ, Q and the mass and energy potentials Mi, Ei are 

Ml = nil cosh£ + 7Ti2 sinh£, 
M2 = Tn2 cosh£ + mi sinh£. 

El = El coshj' +e2sinhj', 

£■2 = £2 coshj'-l-eisinhj'. (A14) 

Eq. (|A12[) depends on standard Pauli-Dirac representation of gamma matrices in block forms (see Eq. (2.28) in [Ij] 
for their explicit forms) and wherq]£| 

Si = 75iA:7±i- (A15) 

2. Vector potentials Af^ in terms of the invariant A{r) 

Comparing E g. (|A12[) with Eq. ([T4)) we find that the spin-dependent electromagnetic-like vector interactions of 
Eq. dill) are [H, HI 

i^^ = {{ei - Ei))P^^ + (1 - GK - ^5G • 7272^ 

At^ = ((£2 - E2))P^ - (1 - G)p>i + '-dG ■ 7i7r, (A16) 



^^ In short, one inserts Eq. l|A3|l into ||A2|I and brings the free Dirac operator ||A4|| to the right of the matrix hyperbolic functions. Using 
commutators and cosh A — sinh A = 1 one arrives at Eq. IIA12II . The structure of these equations are very much the same as that of 
a Dirac equation for each of the two particles, with Mi and Ei playing the roles that m + S and e — yl do in the single particle Dirac 
equation. Over and above the usual kinetic part, the spin-dependent modifications involving GVi and the last set of derivative terms 
are two-body recoil effects essential for the compatibility (consistency) of the two equations 

^^ Just as x^ is a four vector, so is P**. Thus, the time-like and space-like interactions in Eq. | |A7| | become 7^72 and fi ■ 72 only in the 
cm. system due to the fact that from Eq. (IA8I 1. ft = 7? only in the cm. frame. Likewise, S^ = (0, S) only in the cm. frame just as 
is x^ = (0, r) in that frame only. 
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Note that the first portion of the vector potentials is time-like (parallel to P^) while the next two portions are 
space-like (transverse to P^). The spin-dependent scalar potentials Si are 

5i = Ml - mi - ^G72 • dC, 

3-2 = Ah - 7712 + ^G7i ■ d C. (A17) 

We have chosen a parametrization for the vertex invariants C, J^ = —Q that takes advantage of the Todorov effective 
external potential forms and at the same time will display the correct static limit form for the Pauli reduction. The 
logic of the choice for these parametrizations is strengthened by the fact that for classical [2^ or quantum field 



theories [2J] for separate scalar and time-like vector interactions one can show that the spin independent part of the 



quasipotential $ involves the difference of squares of the invariant mass and energy potentials 

Mf = m} + 2mn,S + S^; Ef = e} ~ 2e.^A + A^ , (A18) 

so that 

Mf - Ef = 2m^S + 5^ + 2e.^A - A^ - b^{w). (A19) 

Eqs. P^ and (|A12p involve combined scalar and electromagnetic-like vector interactions (without the separate 
time-like interactions this amounts to working in the Feynman gauge with the simplest relation between space- and 
time- like parts, see Eqs. (jA10[) . (|A11[) . and [l3|, |4J]). In that case the mass and energy potentials in place of Eq. 
(|A18|) are respectively 

Mf = m"^ +exp{2g){2m^S+S^), 
Ef^eM2giA))i{e,-A)'), 
Mf - Ef = exp{'^g{A))[2m.u,S + S^ + 2e^A - A^ - 9{w] (A20) 

so that from Eq. (|A14|) . 



exp(£) = exp(£(^, A)) = ^< + ^^P(2g)(2»^..^+^^) + ^^i + exp(2g)(2777.5+5^) ^^^^^ 

7771 + 7772 

exp(J') = cxp(-C7) 
with 

exp(2^(^)) = (^^^ - G^ (A22) 

or 

-gJ- log(l - 2Alw) - log ^^+^\ (A23) 

2 w 

and the spin-independent minimal coupling appears like 

$5/ = 2777^5 + S'^ + 2ey,A - £■ . (A24) 

3. Interaction terms in the equivalent Relativistic Schrodinger Equation [Eq. (|17p ] 

The Klein-Gordon like potential energy terms appearing in the Pauli form f fTTj) arise from (see Eq. (|A20[) ) 

Mf - Ef = exp(2e) [2777^^5' -h S^^ + 2e.^A - A^ - \?{w)\. (A25) 

To obtain the simple Pauli form of Eq. (fT6|) and the subsequent detailed form in Eq. (|17p involves steps similar to 



those used in the Pauli reduction of the single particle Dirac equation |15j but with the combinations 4'± = tpi ^ 4^4 



and x± = ■02 i 7/^3 instead of the upper and lower components of the single particle wave function. This reduces the 
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Pauli forms to 4 uncoupled 4 component relativistic Schrodinger equations [Tj, [23, |30, 13 • ^^ work in the cm. 
frame in which P = (1, 0) and f = (0,r). We also define four component wave functions il'±,'r]± by [29| 

(l>± = exp(J^ + /C<Ti-f(T2-f)-!/'-|- = (exp J^)(cosh/C + (Ti-fcr2'f sinh/C)'(/'±, 

X± = eyip{F + Kcri-rcr2-f)rj± — (exp J^)(cosh/C + cri'fcr2'f sinhA^)7y±, (A26) 

in which 

1 V 

Z £2"^! + £1^2 

V=E2Mi+EiM2, 

^ = ^^±^. (A27) 

The substitution (|A26[) has the convenient property that in the resultant bound state equation, the coefficients of the 
first order relative momentum terms vanish. 

Using the results in [2^ and [la] we obtain for the general case of unequal masses the relativistic Schrodinger 
equation ( [T7)) that is a detailed cm. form of Eq. ([TB| . In that equations we have introduced the abbreviationJn 

2(J-' + l/r)(cosh2A:-l) ^,2 ^,2 2/C'sinh2/C „. 
r r 

J" (J'' + l/r)(cosh2/C-l) A:'sinh2/C 

'^so = \ , 

r r r 

'^SOD == (l' cosh 2/C - q' sinh 2/C), 
^sox == (<7'cosh2A: + /' sinh 2/C), 

2/C'sinh2A: 2(J"' + l/r)(cosh2/C-l) 2J''/C' V^/C 



$ss == n{r) 



3r 3r 3 3 



$, = ![.(.) + (3-^--/C- + 3A)sinh2/C ^ (^--3 ^ + 1/.) (cosh 2/C - 1) ^ ^_^, ^, _ ^2^j^ 
3 7' r 

^, cosh 2/C -1 /C' , (^' + l/r) sinh 2/C 



where 



n{r) = v^ic - Uj^g + ^^^ ^^^^' + ^'g' - 2^'k;', 

m(r)--iv2g+^g'2-/C'2 + g'^', 

^, ^ (£-g)^£;2M2-^iMi 
2r £;2Mi+£;iM2' 
(£-g/i^a/2-i^2M, 
^" 2r E2M1 + E1M2' ^ ' 

(The prime symbol stands for d/dr, and the explicit forms of the derivatives are given in Eq. (jA30|) ). For L = J 
states, the hyperbolic terms cancel and the spin-orbit difference terms in general produce spin mixing except for equal 
masses or J = 0. For ease of use we have listed below the explicit forms that appear in the above $s in Eqs. (|A28p 
- (|A29p in terms of the general invariant potentials A{r) and S{r). The radial components of Eq. (|T7)) are given in 
Appendix B. 



^* Minor misprints of the equations below have appeared in appendices in [15J and [l4| . The ones presented here are corrected. 
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4. Explicit expressions for terms in the relativistic Schrodinger Equation (|17p from A{r) and S{r) 

Given the functions A{r) and S{r) for the interaction, users of the relativistic Schrodinger equation P7|) will find it 
convenient to have an explicit expression in an order that would be useful for programing the terms in the associated 
equation (jA28p . We use the definitions above given in Eqs. (|A20l ')- (jA22|) . and (|A27p . In order that the terms in Eq. 
(|A28p be reduced to expressions involving just A{r), and S{r) and their derivatives, we list the following formulae 

J,, ^ {C -g'){E2Al2 + EiMi) _ 
2{E2h'h + E1M2) 



W-2A' 



M[_M^_ w (S'{m^ + S) , {2m^S + S'^)A' 



Also needed are 



M2 Ml M1M2 \ W-2A {w-2AY 

X' = ^I^11, (A30) 



uor- I f {ei + e2){Mi + l\l2) ^ {mi + m2){Ei + E2) 
cosh2/C — 



2 V(toi + m2){Ei + E2) (£1 + £2)(Mi + Ah) 
sinh2^ . i (^ (fi+f^KMr + AT^) _ ("^ + m.)(i^r + J^.)', 

2 V (mi + m2)(Si + E2) (£1 + £2)(A/i + Ah) ' ' ^ ^ 



and 



72 



(v^/:-v^g)(i;2M2 + i;iA/i) ^, ,.2 K-mi)^ 

2{E2Ah+EiAl2) 2{E2Ah+EiAh)^ 



^2^^_^'HM? + A4l) 



M1M2 

' V^S-Cto^^ + S) + S"2 45" (m^ + S')A' + {2m^S + S'^)V'^A A{2m^S + S^)A''^ 



M1M2 \ W-2A (u) - 2A)2 {w - 2A)3 

V2g = ^^-^ + 2g'2. (A32) 

The expressions for K(r), m{r), and n(r) that appear in Eqs. f IA28[) ) are given in Eqs. (|A29|) . They can be evaluated 
using the above expressions plus 

V^IC ^ Z!^±Z!£. (A33) 

The only remaining parts of Eq. (jA28[) that need expressing arc those for /' and q' . Using Eq. ()A27P they can be 
obtained in terms of the above formulae. 



Appendix B: Radial Equations 

The following are radial eigenvalue equations [la, [23 corresponding to Eq. p7|) . For a general singlet ^Jj wave 
function u^sj = '^JOJ = '^o coupled to a general triplet ^Jj wave function uju = ui, the wave equation 

(f2 J(J+l) 
{-TT + 9 + 2m«.5 + S^ + 2e.u,A - A^ + $z3-3$ss}uo 



= 6'uo, (Bl) 
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is coupled to 






- 2$so + ^ss + 2$r - 2$5ot}mi + 2v/J(J + 1)($sod + $sox)uo 

= 6^wi. (B2) 

For a general 5"=!, J = -L + l wave function uj^uj = u+ coupled to a general S'=l, J = L — 1 wave function 
uj+iij = U- the equation 



d^ JiJ - 1) 
{ — r^ + ^ — - + 2™«.5' + 5*2 + 2£^„A - A^ + $£, 

2(J-1), 
+ 2( J - l)$so + *ss + ^7^(*soT - $t)}u+ 



2J+1 ^ ^^^ " ^^"^ + 2)$sot}u_ 
6^^+, (B3) 



is coupled to 



^ , (J + l)(J + 2) 



{-TT + ^ ^^ + 2"i«i'S' + 5^ + 2e„A - A^ + $ 



2f J + 2) 
2(J + 2)$so + *ss + -^J^^^soT - $t)}w- 

^^^^{3$^ + 2( J - l)<i>soT}^+ 

:6^t_. (B4) 



For the uncoupled ^Pq states the single equation is 

d2 2 



{--7^ + ^ + 2"i«;5' + S-^ + 2£.„A - A^ + $ 



D 



4$so + *ss + 4($soT - $t)}m_ = 6^M_. (B5) 



1. Specialization to vector interactions, equal masses and J = Q. 

In this case we need only consider the ^Sq and '^Pq states. The corresponding equations are 

{-^ + 2e-^ - ^' + *d-3$ss}mo = &'mo, (B6) 



and 



- 4$so + *ss + 4($50T - *t)}w- = 6^M_. (B7) 
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Wc consider the explicit forms for the quasipotentials given above that appear in these equations for the case of vector 
interactions only, for J — and equal masses. In that case we have 



g' 
a 

J' 

X' 


3^' 
2 ' 


W-2A' 
= 0, 


■W-2A 

{C - J') g' 

2 2 ■ 



Also needed are 



and 



v'^T^-^v^g, 



(B8) 



1 1 

cosh2/C = cosh^ = -( , — + Jl - 2Alw), 

2\/l-2A/w 

sinh2A: = -sinhg = --( , - - Jl~2Alw), (B9) 

2 ^1 - 2A/w 



V^C = 0, 



V^A 



v^j=-v^g = —-2g'^. (BIO) 

w ~ 2A 



In that case we have that the combination for the ^^o equation is 

•^D - 3$ss 

^ 2(J-^ + l/r)(cosh2/C-l) ^ ^,2 ^ ^,2 ^ 2rsinh2/C _ ^^^ ^ ^^^^^ 

2rsinh2/C , 2(^- + l/r)(cosh2/C-l) „ww^„2;^ 

— 3K(r) 1 2y^ /C + V aL 

r r 

= v2(- J- + /c-f - g - /c) + J-" + ^g'2 + -iT'g' 

2 4 

= 0, (Bll) 

while the combination that appears in the '^Pq equation is 

^D - 4$so + ^ss + 4($soT - *t) 

2(J-' + l/r)(cosh2/C-l) ^,2 ^,2 2A:'sinh2A: .^ 

r r 

AT' 4 ( J"' + 1 /r ) (cosh 2A: - 1 ) AJC' sinh 2A: 

r r r 

,, 2A:'sinh2/C 2(J-' + l/r)(cosh2A: - 1) 2J-'/C' V^/C 
+ ,(,) + _ _ + _ _ 

,^,cosh2/C-l 4/C' 4(J"' + l/r)sinh2/C 

_ 4 ^ (3.F'-/C' + 3/r)sinh2/C ^ (^' - 3^ + l/r)(cosh2/C - 1) ^ 2_^,^, _ ^2 

3 r r 

8A' / A' \ ^ 2V^A 



r{w-2A) ' "V'^-S^y ' w-2v4' ^ ^ 
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Thus wc have the two J = single component equations reducing to 



and 



{-j^ + 2e^A~A^}ua^b'uo, (B13) 

dr^ 7-^ 7' (u; — 2A) \w — 2^/ w — 2A 

We consider the case in which 

r 

V^A = 47rJ(r). (B15) 

In that case 

8^' 8a 4 



r (w — 2A) r^ [wr + 2a) r^o r^ ' 

o^ ^' \' 8 / a \' 2 , , 

Vw-2A/ r2 \wr + 2a/ r->o r^ 

This displays explicitly how the spin-orbit and other effects completely overwhelm the angular momentum barrier 
leaving a nonsingular potential at the origin In particular, combining with 2/r^ we obtain 



8A' / A' " ^ 



J.2 r{w-2A) \w-2A) {r + 2a/wY 

From this we obtain Eq. 



(B17) 



2. Specialization to vector interactions, equal masses, and J — L > 0. 



In this case we need only consider the ^Jj and '^Jj states. The corresponding equations are 



i-T^ + '^^'^■^^^ + 2£-A -A^ + $z3-3$ss}uo = 6'uo, (B18) 

and 

{-TT + t + 2e«;A -A^ + ^D- 2$so + '^ss + 2$t - 2$sot}ui 

=: fe^Ul. (B19) 

The first equation simplifies as before $£i=3$ss while for the second equation we have 

^D - 2$so + ^ss + 2*T - 2$sOT 
2» , -^In ni'i 



^^ + V"g - Q' 
r 

2 A' „ / A' \ ^ V'^A 



rw-2A \w-2AJ w - 2A 
Hence, our two J ~ 1 uncoupled equations become 

rf2 2 



(B20) 



{-^ + ^ + 2en,A-A^}uo^b^uo, (B21) 



and 
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= b^ui. (B22) 

Appendix C: Solutions of Eq. (|3ip for Usual and Peculiar ^So Bound States 



Let us use the Coulomb variable r — xJEwOl so that our ^5*0 equation becomes 

u = 2;'^+^w(a;)exp(— fta;), (CI) 

in which the two solutions for A are 

A+-i(-l + v/l-4a2), 

A- = ^(-l-v/l-4a2). (C2) 

corresponding to the usual and peculiar solutions respectively. Then our equation becomes 

- ." + 2./. - ^iA±lK + MA±I> _ 2 , ^ 0^ (C3) 



x x 



Let 



and we obtain 



^;=^7;„^a:"^ (C4) 



nr=0 



For bound states we have 



We let 



(2Kn, - 2 + 2k(A + 1)) ,^_, 

(nr + l){nr + 2(A + Ij) 



^ n,. =0,1,2, .. (C6) 



n,. + A + 1 



?i' = ?i,. + A + L (C7) 

If A were an integer then this would be the principle quantum number n. We write 

i-^---^ + n^)u^O, (C8) 



dx 



2 9 

'^ X X 



d^ 1 a^ 1 
ay yn y^ 4 

where x = y/ (2k) , so that [ig] 

u{y) = exp{-y/2)y'+'Lll+\y) (CIO) 

Let 

^ = ^ = :r^, (Cll) 
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A+l 



and so our radial wave function is 

n' \ n' J "^ n' 

The corresponding hydrogenic radial wave function is 



nao \naQ J nap 



Using the result 146| for the hydrogenic wave function 



0" r„2 



[n^ - 5L{L + 1) + 3] 

u 

and identifying L{L + 1) — > — a^, n ^ n' ,ao — > 1/ {SwOt) we see that for our states 



(0 = 



6 {ewa)' 



;[n'^ + ha^ + Z]. 



(C12) 



(C13) 



(CM) 



(C15) 



Our total cm. energy eigenvalues come from 

(4 - ml) 






£«,(! + ^) = w„, 



-,12 



?l' = ?lr + A + 1, ?lr = 0, 1, , 



(C16) 



In the static limit case for which m2 >> rni we use w = m2 + e in which e << m2 includes the rest mass and 
binding energy of particle 1. Then 

?77iTO2 

rriw = > nil, 

1112 +e 

m\ + 2em2 + £^ — mf — m| 

2771 1TO2 

2em2 + £^ — m? 
-)► 

2toi7TT,2 

In that case the above solution would be for the binding energy 

ra 

e = ± 



(1 + ^) 



Since we do not include negative energies we dispense with the lower sign. 

Let us solve for the total cm. energy in the case of equal masses mi = 7TI2 = rn, 






w"^ - 2to2 

T- = /(a) = 



1 



2to 

2 r,™2 



(1 + S) 



Thus the solutions are 



W± — V 2771 



1 



\ V( 



(n,+A±+l)^ 



(C17) 



(C18) 



(C19) 



(C20) 
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Since i = we take our principle quantum number to be n = n,, + 1. This leads to the results in the text for the 
spectrum. The value of (r^) for the peculiar ground state is 



r')- 



r[n'2 + 5a^ + 3] 



6 (ewCt) 

8 (e^a) 4 



^ ^ ^ ^-^- = ^ (C21) 



so that y/ir"^)- is the electron Compton radius. For all of the usual states and the remaining peculiar states they 
have the following forms 

{r')+ = ^^Kf + 5^2 + 3] = lli±^[(n + X+)^ + 5a^ + 3], n = 1, 2, 3..., 

,2\ _ "- r^/2 , K^2 , oi _ (n + X-) 



r" _ 



5^2+3] = ^ — '. '-^[{n + X_y + 5a^ + 3], n = 2,3..., (C22) 



i-y j ' *■ I ^\-X^ I ^J ' t-\ 



and we see that the size of the nth usual state is very nearly the same as the size of the n + 1st peculiar state. 

In light of this one might wonder how the excited peculiar states (which have the size of angtroms) can be orthogonal 
to the peculiar ground state, that has size of a Compton wave length. As an example, as seen from Eq. ( IC5|) the 
first node of the first excited state occurs at 

a; = (A_ + l)(A_+2)-a^ 

r^^^y^. (C23) 

Ew m 

which is on the order of 560 fermis. 

Appendix D: The Connection bet^veen F\{r],br) and G\(r],br) 

We begin with [13-113 

Fa(p) = CA(r/)p^+^ exp(-zp)il/(A + l-iri,2X + 2; 2ip), (Dl) 

and 



exp(— i7rA/2) 



Ga(p) = 77 |r(A + l + ^^)\ cM^vm[ ^,,\ , : L W,,.x+i/2 i^ip) 



+ 



^^M^-.„„A+i/2(-2*p)]. (D2) 



r(A + l-i77) 
We introduce the Coulomb phase shift 

<yxiv) = T^[log(r(A + 1 + .7;) - log(r(A + 1 - Z77)], 
Zi 

r(A + 1 + ir/) = |r(A + 1 + 177)1 exp(zfTA(ry)), (D3) 

and so 

Ga(/o) = -[exp(-i[tTA(?7) - (A - iT])Tr/2)Wi,,,x+i/2 C^ip) + exp{i[ax{7]) - (A - iv)T^/2)W-i^,x+i/2 {-2ip)] 

= i[7A-(A,r7,p)+V'+(A,r7,p)], (D4) 
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where 

V'+(A, V, P) = exp(i[crA(r?) - (A - Jr/)7r/2])W^_,^,A+i/2 (-2ip), 

'Z'- (A, ??, /o) =exp(-i[crA(?7)- (A-i7;)7r/2])VKi,,,A+i/2(2ip), (D5) 

and since X,r], p are all real 

^.{X,V,p)=r+i\v,p)- (D6) 

Also we have 

^A (p) = ^ [V'+ ( A, 77, p) - V- (A, ry, p)] , (D7) 

V'±(A,77,p) = Ga(p)±iFa(p). 

Note that since the Whittaker function Wk,.v{z) is an even function of /i we have that 

V'+(-A- l,ihP) = exp(i[(T_A-i(?7) - (-A- 1 - ir])Tr/2])W^„^,^x~i/2{-2ip) 

= exp{ix{X,r])) exp{i[axiv) - (A - «?7)7r/2])M^_j^,A+i/2 (-2ip) 

= cxp(ix(A,77))'0+(A,77,p), (D8) 



where 



Similarly 



As a result of this we have 



a;(A, 7/) = (A + ^)tt + a_x-i{v) - '^xiv)- (D9) 



^_ (-A - 1, ry, p) = cxp(-ix(A, 77))^'- (A, ry, p). (DIO) 



1 



F_A_i(p) = -[V^+(-A-l,77,p)-^._(-A-l,r,,p)] 

= TT [exp(ia;(A, 77))?/'+ (A, 77, p) - exp(-ix(A, 7y))'0_ (A, 7/, p)] 

27 

= l[exp(7a;(A,7;)) [Ga(p) + iFA(p)] - cxp(~7x(A,77)) [Gxip) - iF^ip)] 
= cosx{X,r])Fx{p) + smx{X,'n)Gxip), (Dll) 

and thus 

= ^-^-^(^)-^';^^^(^'^)^^(/^) . (D12) 

sinx(A,7]) 

Appendix E: The Variable Phase Method of Calogero 

Here we outline the variable phase method, first applied to short range potentials and then to long range potentials. 
We begin with the short range potentials. We consider the following two sets of differential equations 



u" + {b^ - W)u = 


= 0, 


S'/ + (^2 - Wl)u, -- 


= 0, i = l,2, 


MO) = 


= u{0) = 0, 


Wi = 


L{L + 1) 
^2 ' 



(El) 

where W{r) is a short range potential less singular at the origin than const /r^ and 

^1 (^) ~ ]L (br) — > const sin(6r — Ltt /2) , 

722(7') = —nL{br) -^ const cos{br — Lit/ 2). (E2) 
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Let 

u{r) = a{r){cos S L{r)ui{r) + am5L{r)u2{r)) 
u[r — > oo) = cortsf(cos(5i(r — > oo) sin(6r — i7r/2) 
+ sin(5i(7- — > oo) cos(for — L'k/2) 

= constsin{br - Lt:/2 + 5l{oo)) -^ sin(6r - L-k/2 + 5l), (E3) 

and so 

5l = 5l{oo). (E4) 

To find the differential equation that (5l('') satisfies, define 

u {r) = a(r)(co8 5L(r)ui{r) + SYii5L{r)u2(r)), (E5) 



and so 



Then 



u'{r) {cos 5 L{r)u'i{r) +SUY 5 L{r)u'2{r)) {u'i{r) + tan 5 L{r)u'2{r)) 

u{r) {cosdL{r)ui{r) +smdLir)u2{r)) {ui{r) + ta.n6L{r)u2{r))' 

tan 6L[r) = ._..,.. ., , . , .. ■ (E6) 

[U2[r)u[r) — U2[r)u[r)) 



S'^ir) sec^ SUr) = <5l(0(1 + tan^ dUr)) 

{U2U' — U2U) (u'lu — Uiu") — {u'lU — Uiu') {u2u" — u'J^U) 



(u2(r)u'(r') — U2(r)u(r))^ 
{W - Wi)u% 



and so 

Further manipulations lead to 



{W — W7)QL(6r)eos(5L(r) — hL{br) sin(5L(r))^ 



Note that in case of type two reference potentials {W = W7/(r) = 0) we would obtain 



and with 



we obtain 



(E7) 



{U2U' — U^u)"^ ' 

where we have used the Wronskian relation 

U2u'i — UgUi = const = 6, (E8) 



{W~Wi)b 

odn^ TTT — 77 ^TTI- (^9) 

sec^ 0l{U2u' /u — U2) 



(ElO) 



tan7LM - ._^) )' 1}.\ . , (Ell) 

{U2{r)u (r) — U2{r)u{r)) 



ui{r ^- 0) ^ br, 
U2{r^0)^l 
u(r ^ 0) = c(6r)^+\ (E12) 



^^,br)^+^ ^ brcb(L + l)(br)^ 
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and so wc obtain the same boundary condition as with the type I reference potentials. From Eq. (jlOSp 

7iW=-ysin2(fer + 7L(r)), (E14) 



at short distances becomes 

7l(0) = - ^^^ ' ^^ sin^(& + 7^(0)), (E15) 



L{L+\) . 2^ 



with the solution given in Eq. (|109|) . 

Next we sketch an analogous derivation for the phase shift equation which involves long range potentials corre- 
sponding to Eq. (j88p in which the Coulomb potential appears. As discussed in the text we begin with the following 
two sets of differential equations 

u" + {h^ - W)u = 0, 
u'l + {b'^-Wui)u,=Q, i = l,2, 

2£tua a^ 



Win 



r 



W= , ' ,. -^^-4- (E16) 

Note that the total potential plus barrier term W appears in the equation for u. We are not including the angular 
momentum barrier in the definitions of Ui (r) . 
The solutions ui,'S2 to 

+ (^2 ^ 2£^ ^ ^)-^ = 0, * = 1, 2, (E17) 

are Coulomb wave functions 

ill = aF\^ + cG\^ 

U2 = dFx^ + fGx±. (E18) 

We choose the constants so that ui has the same behavior at the origin that u does. 

Even though four functions are listed here, only two are linearly independent (see Eq. (|80|). To determine the 
phase shift equation let us write down first the wave function u{r) in terms of ui,U2 

u(r) — a{r)(cos^±(r)ui{r) + sm^±(r)u2{r)). (E19) 

In that case 

u(r — > oo) — > (cos7±(r — ?► oo) sin(&r — 77log26r + cta± -^ A±7r/2) 
-|-sin7±(r — ?► oo) cos(6r — 7]log26r + ctaj- — A±7r/2) 

= sin(fer - 77log2for + ax± - A±7r/2 + 7±(oo)). (E20) 

This defines the phase shift function 7±(?') and its relation to the asymptotic behavior of u{r). On the other hand 
since u(r) is the wave function for a potential that includes at r >> 2a /w the modified angular momentum barrier 
(2 — a^)/r^ = k(k + l)/r^ in addition to the Coulomb term, we must have 

uir -> oo) -> sin(6r - ?7log26r + a,^- k,'k/2 + d'^), (E21) 

and so comparison gives 

o-A± - A±7r/2 + 7±(oo) =(j^- ki:/2 + 5^,. (E22) 

Thus with 

K(K + l) = 2-a^ (E23) 

-l + \/9-4a2 
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the full phase shift is 

Sk + o-k = o-a± +{k- A±)7r/2 + 7±(oo) 

= argr(A± + 1 + ir/) + (k - A±)7r/2 + 7±(oo) 

~ argr(A± + 1 + «/) + (1 - A±)7r/2 + 7±(oo). (E24) 

To find the differential equation that 'y±(r) satisfies, define 

u'{r) = a{r){cosj±(r)u'i{r) + sin7±(r)M2(r)), 
u{r) = a{r){cosj±(r)ui{r) + smj±{r)u2{r)). (E25) 

Then following a procedure similar that given in Eqs. (jE5[) we obtain 

, , u'-,(r)u(r) — ui(r)u'(r) /^„„x 

tan 7± (r) - .-^ /. _,. \ , ■ (E26 

(U2[r)u'[r) — U2[r)U[r)) 



Also 



7i (0 sec^ 7± {r) = - ,_ ^"'_^, „ , (E27) 



where we have used the Wronskian relation 

U2u'i — u^ui — const 

= cos() cos()(5 - -) + sin() sin()(6 - -) 
r r 



and so 



{W - Wni)h 
sec^ 7±(u2m'/u — Uj)^ 



Now use 

u'(r) (cos7±(r)'Uj(r) + sin7±(r)M2(r)) 
u(r) (cos7±(r)Mi + sin7±(r)-U2) 

and hence, with ui = F\j_ , U2 = G\j_ we have 



(E28) 



7^(0 = - 2^;: f-M2 - (E29) 



(E30) 



,,, (^-V^^//j)(cos7±(0J^A±(r)+sin7±(r)GA^(r))^ 
7±(r) = ^ . (E31) 

Because of the 2/r^ behavior of W for large r one will have to integrate quite far to obtain a convergence for 7± (r) 
and after that one must subtract the phase shift — 7r/2 due to the 2/r^ angular momentum barrier. An alternative 
form of this equation is 

tan 7±(r) = . (E32) 

The question now arises about the boundary condition at the origin for 7± (r) . We focus on Eq. (jE26p to determine 
the boundary condition at the origin for 7±(r), 

, ^ , X u'^{r)u{r) ~ ui{r)u' [r) 

tan7± r) = ,_,,,,, _, , , ,- ■ (E33) 

{U2{r)w (r) — U2{r)u[r)) 

We determine the behavior at the origin by evaluating the right hand side for very small r. The dominant term for 
the quasipotential for both case is —a^ /r^ . Thus it is sufficient to focus on the first case 
We use Eq. (p33)) with 

— ^ u = 6^u, (E34) 



dr"^ (r + 2a /w)^ r r 
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and 

d'^ 2e.ii,a a^ 1 _ ,9- 

TT 7 r "1,2 = ""1,2, 

a?'^ r r^ J 

M2=GA(&r)- (E35) 

At short distance, the potential energy for u is the same as that for mi_2- At very short distance, we choose 

u, ui — > const — > r ^ , 
M2 -^ ar^^^ +br^^. 

Clearly then tan 7± (0) = as the numerator vanishes in both cases where as the denominator is proportional to the 
Wronskian of ui and U2 which is b^. This case allows us to integrate either Eq. (jE31[) or (|E32p with the boundary 
condition of tan7±(0) = 7±(0) = 0. 

To find the total wave function we need to find the additional differential equation for the amplitude of the wave 
function. We use 

u'{r) — a{r){cosj±{r)u'^{r) + sin7±(r)'U2(r)) 
= a'(r)(cos7±(r)-ui(r) + sin7±(r)M2(?')) 
+ a(r)(cos7±(r)u']^(r) + sin7±(r)u2(^)) 
+ j'^a{r){—smj±{r)ui{r) + cos7±(r)u2(r)), (E36) 



and thus, using Eq. (jE31|) 

a'jr) ^ Wir)-Wiii {ul-ul)sm2j±{r) 
a{r) b ^ 2 



U1U2 cos27±] 



i , f ^ } I WJr) - Win jr) , {ul - ul) sin 27± (r) _ _ ^ . ... ,p„„. 

a(r) = a(ro)exp{- / — ^^ — ^^[^ '— U1W2 cos27±(r)]}. (E37) 



So, the total wave function is 



M(r) = a(ro)exp{- / [^ '— M1U2 cos27±(r)]} 

X (cos7±(r)ui(r) + sin7±(r)fZ2(?'))- (E38) 
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